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1 Introduction 



Our aim in this article is to describe some links which are starting to emerge between two 
previously-separated areas of mathematical physics - the theory of integrable models in 
two dimensions, and the spectral analysis of ordinary differential equations. 

The study of integrable lattice models has been an intriguing part of mathematical 
physics since Onsager's solution of the two-dimensional Ising model. Lieb and Suther- 
land's work on the six-vertex model showed that this was not an isolated phenomenon, 
and with Baxter's work at the beginning of the 1970s the full richness of the field began to 
be appreciated by a wider community. Since that time interest in the subject has grown 
steadily, receiving a particular boost of late from the links which exist with integrable 
quantum field theories. Many different methods exist for the solution of these models, 
and a technique which will be very important in the following goes by the name of the 
'functional relations' approach. The idea, initially put forward by Baxter, is to show that 
quantities of interest satisfy functional relations. When combined with suitable analytic- 
ity properties, these relations can be highly restrictive and often lead to exact formulae 
for quantities of physical interest. 

In a parallel chain of development which also dates back at least to the early 1970s, 
Sibuya, Voros and others have shown that functional relations have an important role 
to play in a rather more classical area of mathematics, namely the theory of Stokes 
multipliers and spectral determinants for ordinary differential equations in the complex 
domain. However, it is only recently that the existence of a precise link - an 'ODE/IM 
correspondence' - between these two areas has been realised. The aim of these notes is to 
provide an elementary introduction to this connection and its background. For most of 
the time the focus will be on the simplest example, which connects second-order ordinary 
differential equations to integrable models associated with the Lie algebra SU{2). 

We begin, in the next section, with a short introduction to the types of spectral 
problems which will be relevant, mentioning in the process our third main theme, the 
intriguing reality properties of certain non-Hermitian spectral problems which arise in 
the study of "PT-symmetric' quantum mechanics. Integrable lattice models and their 
treatment via functional relations are introduced in section [Sj we also discuss briefly the 
recent development of these ideas within quantum field theory. The differential equations 
side of the story is explained in section U after which the link with integrable models is 
made precise in section [H Some applications and generalisations of the correspondence 
are outlined in section [U and section [7] contains our conclusions. Various pieces of 
background material have been collected in the appendices, including an explanation of 
the algebraic Bethe ansatz in appendix Rl 

More on the ODE/IM correspondence can be found in references [1-19]. All of this 
work rests on earlier studies by, among others, Sibuya [20], Voros [21], and Bender et 
al [22-24] (on the ODE side) and by Baxter [25,26], Kliimper, Pearce and collabora- 
tors [27, 28] , Fendley et al [29] , and Bazhanov, Lukyanov and Zamolodchikov [30-32] on 
the integrable models side. 

We have aimed to make this article accessible to readers with backgrounds in both 
the integrable models and the differential equations communities; for this reason, we have 
tried to keep the treatment relatively elementary. More details can be found in many 
places [20,21,33,34]. 

Those readers primarily interested in ODEs may prefer to concentrate on section [21 
briefly read sections 13.21 13.91 and 13.101 to get a flavour of the integrable model picture, 
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then move on to sections HI [S] and El while those primarily interested in integrable models 
may like to look initially at sections [21 [3l H] and then at 15.11 15.21 and 15. 6i 
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2 Prelude: three reality conjectures in PT-symmetric quan- 
tum mechanics 



Spectral problems, some of a rather unconventional nature, will play a central role on the 
'differential equations' side of our story. Rather than launch straight into technicalities, 
we shall warm up in this preliminary section by describing an intriguing class of problems 
much studied by Bender and collaborators and others in recent years. It all begin^ with 
a question posed by Bessis and Zinn-Justin, sometime near 1992: 
Question 1: What does the spectrum of the Hamiltonian 

n = p^ + ix^ (2.1) 

look like? 

This is a cubic oscillator, with purely imaginary coupling i. (Strictly speaking, Bessis and 
Zinn-Justin, motivated by considerations of the Yang-Lee edge singularity [42-44], were 
initially interested in more general Hamiltonians of the form + + igx^, from which 
the above problem emerged as a strong-coupling limit.) The corresponding Schrodinger 
equation is 

- T^i^ix) + ix^Mx) = Eipix) (2.2) 

and we shall declare that the (possibly complex) number E is in the spectrum if and only 
if, for that value of the equation has a solution ^ipi^) on the real axis which decays 
both at x ^ — oo and at x — > +oo, as illustrated in figure [Tpl . Notice that the differential 




Figure 1: A wavefunction decaying at x = ±00. 

equation forces the wavefunction ip{x) to be complex, even for real values of x and E. 
And since the Hamiltonian is not (at least in any obvious way) Hermitian, the usual 
arguments to show that all of the eigenvalues must be real do not apply. Nevertheless, 
perturbative and numerical studies led Bessis and Zinn-Justin to the following claim: 

Conjecture 1 [47]: The spectrum of W is real, and positive. 

What might be behind this strange property? Bender and Boettcher [23] stressed the 
relevance of VT symmetry. To be more precise, 'P', or parity, acts by sending x to —x 
and p to —p while T, time reversal, sends x to x, p to —p and i to —i. Note that both 

^While this question initiated the line of work we want to describe here, similar curiosities had in fact 
been observed before - see, for example, [35-39]. For further historical discussions, see [40,41]. 

^To be more precise, the decay should be fast enough that ip{x) lies in L^(R), the space of square- 
integrable functions on the real axis. This means that we are actually discussing the so-called point 
spectrum oi H - see, for example, [45,46]. 
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V and T preserve the canonical commutation relation = i of quantum mechanics 
even if x and p are complex [24] . 

As shown in [24], VT invariance implies that eigenvalues are either real, or appear in 
complex-conjugate pairs, much like the roots of a real polynomial. But, just as the typical 
real polynomial has many complex root^, on its own VT invariance of the Hamiltonian 
does not guarantee reality. This is elegantly shown by the following generalisation of the 
Bessis-Zinn- Justin problem, proposed by Bender and Boettcher [23]: 

Question 2: What is the spectrum of 

nM=p'^ - {ixf^ (M real, > 0) ? (2.3) 



Later, it will turn out that the passage from question 1 to question 2 corresponds to a 
change in a parameter in a lattice model, or equivalently to a change of a quantum group 
deformation parameter in a Bethe ansatz system. But for now, the generalisation is 
appealing because it unites into a single family of eigenvalue problems both the M = 3/2 
case, for which we have the Bessis-Zinn-Justin conjecture, and the much more easily- 
understood M = 1 case, the simple harmonic oscillator. Furthermore, for all M the 
problem is 'PT-symmetric. The Schrodinger equation is now 

- - (ixf^'i^ix) = EMx) (2.4) 

and again we ask for those values of E at which there is a solution along the real x-axis 
which decays at both plus and minus infinity. Two details need extra care: for non- 
integer values of 2M, the 'potential' —{ix)"^^^ is not single- valued; and when M hits 2, 
the naive definition of the eigenvalue problem runs into difficulties. The first problem 
is easily cured by adding a branch cut running up the positive imaginary x-axis. The 
second is more subtle, and its resolution more interesting; it will be discussed in greater 
detail in section S] below. 

For the moment, we shall agree to keep M below 2. Even so, there is an interesting 
surprise. Figure [2] is taken from [4] - it reproduces the results of [23] . Clearly, something 
strange occurs as M decreases below 1. Infinitely-many real eigenvalues pair off and 
become complex, and only finitely-many remain real. By the time M has reached 0.75, 
all but three have become complex, and as M tends to 0.5 the last real eigenvalue diverges 
to infinity. In fact, at M=0.5 the problem has no eigenvalues at all, as can be seen by 
shifting ix to ix — E and solving the resulting equation using an Airy function. For M > 
1, numerical results combined with various pieces of analytical evidence indicated that 
the spectrum was entirely real, and positive, and so Bender and Boettcher generalised 
conjecture 1 to 

Conjecture 2 [23]: The spectrum of Ti.M is real and positive for M > 1. 

The 'phase transition' to infinitely-many complex eigenvalues at M = 1 can be interpreted 
as a spontaneous breaking of VT symmetry [23] . 

One further generalisation of the Bessis-Zinn-Justin conjecture will be relevant later. 
Consider 

famous result of M. Kac shows that the expected fraction of real zeros of a real polynomial of degree 
n with random (normally-distributed) coefficients tends to zero as n ^ oo, as ■|log(n)/n. See [48,49], 
and, for example, [50,51]. 
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Figure 2: Tiu = ~ {ix)'^^ : real eigenvalues as a function of M. 

Question 3: What is the spectrum of 

nM,i = / - {ixf^^ + l{l+l)/x'^ (M and / real, M > 0) ? (2.5) 

This amounts to studying the effect of an angular-momentum-like term l{l+l)x~'^ on the 
Bender-Boettcher problem, and it was first investigated in [4]. Note that we continue to 
impose boundary conditions at x = ±oo, in the way stated just after equation (j2.2p above. 
With the angular-momentum term included we need to specify how the wavefunction 
should be continued around the singularity at x = 0; given the choice to place a branch 
cut on the positive-imaginary x axis this continuation should be done through the lower 
half plane. (There will be much more discussion of boundary conditions later, so we 
won't go into this detail any more for now.) Again, a combination of numerical and 
analytical work gave strong evidence for 

Conjecture 3 [4]: The spectrum of TCmi is real and positive for M > 1 and |2/+1| < 
M+1. 

Although a small angular-momentum-like term does not have a significant effect while 
M > 1 and the eigenvalues all remain real, for M < 1 it can make a remarkable difference 
to the way in which they become complex. Figure [3] shows the spectral plot for I = 
—0.025, and reveals a dramatic change from the earlier / = plot: the connectivity of 
the real eigenvalues has been reversed, so that while for I = the first and second excited 
states pair off, for / = —0.025 the first excited state is instead paired with the ground 
state, and so on up the spectrum. With this in mind, it may be hard to see how it is 
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Figure 3: TLm,i = — (ix)"^^'^ + l{l+l)x~'^ : real eigenvalues as a function of M for 
I = -0.025, /(/+!) = -0.024735. 

possible to pass between the sets of spectra depicted in figures [3] and [2] simply by varying 
the continuous parameter I from —0.025 to zero. The mechanism should become clear 
by looking at the sample of intermediate pictures of figure HI 

A final generalisation allows for an even richer phenomenology. Adding an inhomo- 
geneous term —a{ix)^^~^ to the potential for 7iM,i gives a three-parameter family of 
problems: 

nM,a,i =P^- {ixf^ - a{ix)^-' + l{l+l)/x^. (2.6) 

Again, the first question to ask is whether the spectrum of 'HM,a,i is entirely real. Some 
general results will be described later in this review, but for now we illustrate the situation 
by giving some 'experimental' data for the case M = 3. Special features of this particular 
case make it desirable to trade the parameter I for p := -v/3(2^+l); this being understood, 
figure [5] shows, below the cusped line, the region of the (a,/3) plane where all eigenvalues 
of 7i^^ct,i are real. One last comment is worth making here: the alert reader might protest 
that substituting Af = 3 into p.6p results in a Hamiltonian which is manifestly real, with 
a potential bounded from below, which should surely have an entirely real spectrum for 
all values of a and I (or p). The question is a fair one, but it fails to take into account 
the fact that we have implicitly imposed boundary conditions which analytically continue 
those at M = 1. This continuation takes M past the value 2 at which we already observed 
that there would be subtleties in defining the spectrum. We shall return to this point 
later. 
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Figure 4: Real eigenvalues of — (ix)"^^ + l{l+l)/x'^ as functions of M, for various 
values of /. 

Bender and Boettcher's observation [23] has sparked a great deal of interest in reality 
properties in non-Hermitian quantum mechanics; a (small) sample of related work on the 
reality issue is provided by references [52-112]. In this already-long article, we won't have 
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Fi gure 5: Tl^^a,i '■ region of unreality in the plane, where p — \/3(2/+l). 

space for any further discussion of non-Hermitian quantum mechanics in general; more on 
current issues in the field can be found in, for example, [41], the review article [113], the 
conference proceedings [114], and references therein. Instead we just remark that reality 
properties in PT-symmetric quantum mechanics of the sort described above have turned 
out to be surprisingly hard to establish by conventional means. An interesting byproduct 
of the ODE/IM correspondence has been a relatively elementary proof of conjectures 1, 
2 and 3, and an understanding of many features of figure [5j We shall give this proof in 
section 16.21 below; it relies heavily on certain functional relations which had first made 
their appearance in a very different context: the theory of integrable lattice models. A 
rapid introduction to the background to this material is our next subject. 
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3 Integrable models and functional relations 



In this section we shall introduce the integrable lattice models and quantum field theories 
that will be relevant later. One particular technique for their solution, the 'functional 
relations' approach, will be highlighted. We start with the lattice models, and then 
discuss what is known as the 'continuum limit' in preparation for the link with quantum- 
mechanical problems. 

3.1 Generalities 

Lattice models provide a way to understand the behaviours of magnets, and other sub- 
stances, which exhibit a number of distinct 'phases' depending on the values taken by 
external parameters, such as the temperature. The simplest example is the Ising model, 
where a macroscopic magnet is modelled by a square lattice of microscopic 'atoms', each 
of which can exist in one of two states of magnetisation, up or down. By coupling nearest- 
neighbour atoms by an interaction which favours the lining-up of their magnetisations, 
a system is produced which captures the tendency of real magnets to exhibit an overall, 
macroscopic, magnetisation at low temperatures, which is then lost as the temperature 
is increased past some critical value. At this point, the model is said to undergo a 'phase 
transition', from an 'ordered' to a 'disordered' phase. The qualitative features of this 
behaviour can often be deduced from general arguments - the Peierls argument [115] is 
a good example - but for certain two-dimensional cases it turns out that much more can 
be said, and a number of important quantities can be calculated exactly as functions 
of the external parameters. Crudely speaking, these are the integrable lattice models, 
and the first example found was the two-dimensional Ising model, solved by Onsager in 
1944 [116]. 

Much more can be said on this topic, but for the purposes of this review it is best 
now to move forward by two decades, and to describe the model which will be directly 
relevant to our subsequent story. 

3.2 The six-vertex model and Bethe ansatz equations 

We shall be interested in one particular set of integrable lattice models, called the 'ice- 
type' models, or, a little more prosaically, the six-vertex models. They were first solved 
in 1967, by Lieb [117] and Sutherland [118], and are among the simplest generalisations 
of the Ising model. Good places to look for further details are the book [33] by Baxter, 
and the short review [119] by McCoy. 

The definition of the model begins with an x A^' lattice, with periodic boundary 
conditions in both directions and, in order to avoid some annoying signs later on, N/2 
ever0. Ultimately, we shall take a limit in which N ^ cc. On each horizontal or vertical 
link of the lattice, we place a spin 1 or 2, conveniently depicted by an arrow pointing 
either right or left (for the horizontal links) or up or down (for the vertical links), as in 
figure m 

In the 'ice' picture, each vertex represents an oxygen atom, and each link a bond be- 
tween neighbouring atoms. Sitting on each bond is a hydrogen ion, lying near to one or 
the other end of the bond, according to the direction of the arrow. The 'ice rule' [121, 122] 
states that each oxygen atom should have exactly two hydrogen ions close to it, and two 

^Sometimes the parity of A'^ can be significant, however: see, for example, [120]. 
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Figure 6: A typical configuration of spins for the six-vertex model. 

far away. Translated into the arrows, this implies that only those configurations that 
preserve the flux of arrows through each vertex are permitted. This means that there are 
six options for the spins around each site, or vertex, of the lattice (hence the alternative 
name 'six- vertex'). For real ice, all allowed configurations are equally likely; but if we 
want to generalise slightly, then we can allow for differing probabilities for the different 
options. These overall probabilities are calculated in two steps. First, a number W, 
called a (local) Boltzmann weight, is assigned to each local possibility. If we further 
restrict ourselves to what is called the zero-field case, then the Boltzmann weights should 
be invariant under the simultaneous reversal of all arrows, and just three independent 
quantities need to be specified: 
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= w 
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= a 
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= w 


T 


= b 
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= w 


i 


= c 



(3.1) 
(3.2) 
(3.3) 



The relative probability of finding any given configuration is simply the product of 
the Boltzmann weights at the individual vertices. A first quantity to be calculated is the 
sum of these numbers over all possible configurations - the partition function, Z: 

{cr} sites 

One of the special features of integrable models is that quantities such as the partition 
function (or even better, the free energy per site, defined in equation (j3.9p below) can be 
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evaluated exactly, at least in the 'thermodynamic limit', where and N' both tend to 
infinity. The model under discussion turns out to be integrable for all values of a, b and c. 
Their overall normalisation factors out trivially from all quantities, and it is convenient 
to parameterise the remaining two degrees of freedom using a pair of variables v and r], 
called the spectral parameter and the anisotropy: 

a(z/, r/) = sin{r] + ih') , h{v,ri) = s\n{rj — iv) , c{v,rj) = sin(2?7) . (3-5) 

In calculations the anisotropy is usually held fixed, but it is useful to treat models with 
different values of the spectral parameter u simultaneously. The weights can then be 
drawn as in figure [71 
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Figure 7: The local Boltzmann weights. 



For those familiar with integrable quantum field theories, this picture might suggest 
a relationship between Boltzmann weights for integrable lattice models and S-matrix 
elements for integrable quantum field theories, with the spectral parameter of the Boltz- 
mann weight proportional to the relative rapidity of two particles in the quantum field 
theorjQ. We won't explore this aspect much further here, but a nice discussion can be 
found in [123]. 

One popular tactic for the computation of the partition function employs the so-called 
transfer matrix, T. Introduce multi-indices a = (ai, a2 ■ ■ ■ un) and a' = (a'^, 02 
and set 



, a 



N) 



{ft} 



w 



'u)W 



a2 



'u)W 



03 



{iy)...W 



(3.6) 



A pictorial representation of T(i/) is given in figure [8j The definition involves a sum over 



T(v) = 




Figure 8: The transfer matrix. 

one set of horizontal links, and from the picture it is clear that the matrix indices of T 

^ In fact, in defining the local weights a, b and c we have shifted the spectral parameter by irj from 
the value that would be appropriate for an S-inatrix, so as to give later equations a more standard form. 
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correspond to the spin variables sitting on the vertical links. These can now be summed by 
matrix multiplication, with a final trace implementing the periodic boundary conditions 
in the vertical direction. Thus: 



Trace 



(3.7) 



Calculations can now continue via a diagonalisation of T. Suppose that the first few 
eigenvalues, to > ti > . . . are known, with eigenvectors . . . : 

^TS'*H'^ = t,*k). (3.8) 

a' 

Then, for example, the free energy per site in the limit N' ^ oo can be obtained as 

/ = -^logZ = log Trace [t^'] ~ -^logto • (3.9) 

The eigenvalues to, ti . . . are functions of i' and ??, and the remaining task is to find them. 
This appears to be a very tough problem - T is a 2'^ x 2^ matrix, and quickly becomes 
too large for even the most powerful computers to handle. It is necessary to exploit some 
of the special features of the model, and a popular technique for doing this goes by the 
name of the Bethe ansatz. There are two steps: 

(i) Make a (well-informed) guess for a form for an eigenvector of T, depending on a finite 
number n of parameters i/i, . . . i/„ (the boots'). 

(ii) Discover that this guess only works if the {vi} together solve a certain set of coupled 
equations (the 'Bethe ansatz equations'). 

Letting n vary over a finite range, and for each n taking the finite set of solutions to 
the corresponding Bethe ansatz equations, should then give totality of the eigenvectors 
of T, or at least all those needed to capture the — > oo limit of the systemic. 

The justification of this procedure is an interesting story, or collection of stories, in 
its own right; in appendix A we outline one of the more elegant approaches, a technique 
called the algebraic Bethe ansatz. 

For the six-vertex model, when the dust has settled the Bethe ansatz equations for 
the roots {i/i . . . f„} are 

This is a set of n equations for n unknowns. There is no unique solution, but rather 
a discrete set. For each solution, an eigenvector 1^*) of T can be constructed, with 
eigenvalue 

n n 

t{u) = a^(i/, ri) n gi'^j -1^) + b^i'^, ri) J] ^(^ " '^j) , (3-11) 
3=1 i=i 
where g{iy') := a(i^ — ir],rj)/h{v — ir],r]) = — sin(2r/ -|- iu)/ sin(iz^) . To single out a given 
eigenvector and eigenvalue, supplementary conditions must be imposed on the roots. In 
particular, and this will be important later, the ground state eigenvalue to^i') turns out 
to correspond, in the parameter region < rj < tt/2, Ke(z^) = 0, —r] < Qm{i') < rj, to the 
Bethe ansatz solution with n = N/2 distinct real roots, packed as closely as possible and 
symmetrically placed about the origin (see for example [33,126,127]). This is depicted 
in figure [H 



^We won't go into the interesting question of the completeness of the set of BAE solutions here; 
see [124, 125] for recent discussions. 
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Figure 9: The z^j's are real for the ground state. 



3.3 Adding a twist 

The periodic boundary conditions used above to define the transfer matrix can be mod- 
ified in such a way that integrabihty is not spoilt (see for example [28,128,129]). This 
does not change the free energy per site in the thermodynamic limit, but it does modify 
some subleading effects. 

The twist is introduced by modifying the local Boltzmann weights on one column, or 
seam, of the lattice, say the A^**^ (see figure [TO]l . In the transfer matrix formulation the 
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Figure 10: The twist is introduced by modifying the local Boltzmann weights on one 
column. 



modification amounts to making the substitutions 



W 



a 



'N 



N 



and 
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a. 



a 



N 



N 



a 



'N 



N 



(^) 



(3.12) 
(3.13) 
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in the initial definition (|3.6|) of T. The algebraic Bethe ansatz works almost unchanged, 
with the result that the more general transfer matrix T{i/, (j)) has eigenvalues given by 



n n 



t{u, 4>) = e-*^a^(z., r/) J] g{y^ - v) + e'H'^iu, 7?) [] g{u - uj) , (3.14) 
where the set of roots {z^i . . . satisfy the modified Bethe ansatz equations 

(-irn "t!o"''''*t'''! =-^"'*w^ ■ 

j-J-^ smh(2zry - uj + Vk ) 6^^(ffc, r/j 

3.4 The XXZ model 

There is a well-known connection between classical two dimensional lattice models and 
quantum spin chains. The six-vertex model is related to the (spin 1/2) XXZ spin chain, 
a one-dimensional system of lattice sites with a spin variable taking the value 1 or 2 
at each site, with each spin interacting only with its neighbours. The Hamiltonian is 

1 ^ 

Hxxz = -^Yl {^>J+i + ^M+i - cos 2r? a|c7|+i) , (3.16) 

i=i 

where a" represents a Pauli matrix, 

acting on the spin on the j^^ lattice site. This model is sometimes also referred to as 
the Heisenberg-Ising chain, or the spin 1/2 anisotropic Heisenberg chain. The six-vertex 
twist can be implemented by imposing twisted boundary conditions of the form 

^h+i=^t , af,^,±iay,^, = e^^'H^f±i^!)- (3-18) 

Prom these definitions it is not obvious that the six- vertex and XXZ models should 
be related. However it was noticed in the first work on the six-vertex model [117, 118] 
that its transfer matrix eigenvectors coincided with those of Hxxz , previously studied 
in great detail by Yang and Yang [130] . The initial identification rested on a coincidence 
of Bethe ansatz equations, and was given a more direct explanation when Baxter [131] 
showed that the six-vertex transfer matrix T and the Hamiltonian of the spin chain are 
directly connected through the relation 



d 

Hxxz = -i sm 2r/ — In T{i/] 
dv 



--cos2r/I®^. (3.19) 

-irj 2 



Consequently, if we are able to determine the eigenvalues to^ti . . . of the transfer matrix 
we gain, for free, information on the spectrum of the XXZ model. 

We shall return to a description of the six-vertex spectrum motivated by this connec- 
tion with the XXZ model shortly, but for now we note that for all values of the twist pa- 
rameter, the Hamiltonian (j3.16p commutes with the total spin operator = X^^Li '^f/2- 
Therefore the spectrum splits into disjoint sectors labelled by the spin m = 0, 1, 2 . . . ; 
and the true ground state lies in the m = sector. In the six-vertex model the XXZ 
spin sectors correspond to taking the number of Bethe roots n different from the ground 
state value of n = N/2. The relation between the number of Bethe roots and the XXZ 
spin is n = N/2 — m. 
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3.5 Baxter's TQ relation 

So much, for now, for the Bethe ansatz. There is a particularly neat reformulation of the 
final result, discovered by Baxter, that leads to an alternative way to solve the model. 
The first ingredient is the fact that the transfer matrices at different values of the spectral 
parameter v commute: 

[T(z^),T(z^')] = 0- (3-20) 

(The (standard) proof of this fact is given in appendix A.) Therefore, the transfer matrices 
T(z/) can be simultaneously diagonalised, with eigenvectors which are independent of v. 
This allows us to focus on the individual eigenvalues to(^)) ^\{'^\ ... as functions of v. 
From the explicit form of the Boltzmann weights and the claim that the eigenvectors are 
z^-independent, these functions are entire, and ivr-periodic. 

The second ingredient is the claim that, for each eigenvalue function t(y\ there exists 
an auxiliary function q{y^^ also entire and (at least for the ground state) zvr-periodic, 
such that 

t{v)q{y) = e~''i'a^(y, i])q{iy + 2irj) + eV{u, r])q{v - 2ir]) . (3.21) 

We shall call this the TQ relation, though this phrase should really be reserved for 
the corresponding matricial equation, involving T(z^) and another matrix Q(i^), from 
which the above can be extracted when acting on eigenvectors. At first sight, it is 
not clear why this should encode the whole elaborate structure of the Bethe ansatz 
equations - instead of one unknown function i(i^), we now have two, and all we know 
about them is that they enjoy the curious relationship given by (|3.2ip . But in fact this 
equation, combined with the simultaneous entirety of t{i') and q{i^), imposes constraints 
on t(z^) which are so strong that there is no need to impose the BAE as a supplementary 
set of conditions. (Going further, Baxter was able to establish the TQ relation by an 
independent argument, thereby finding an alternative treatment of the six- vertex model 
which avoided the explicit construction of eigenvectors. He then generalised this approach 
to the previously unsolved eight-vertex model, but we shall not elaborate this aspect any 
further here.) 

The BAE are extracted from (j3.2ip as follows. Suppose, anticipating the final result 
in the notation, that the zeros of q{i') are at vi, . . .Vn ■ Given that g(z^) is ivr-periodic, it 
can - up to an irrelevant overall normalisation - be written as a product over these zeros 
as 

n 

q{u) = JJsinh(i/ - ui) . (3.22) 

1=1 

From (|3.2ip . ti^v) is fixed by q{i^), and from (13.22p . q{i') is fixed by the set {fj}. To 
determine the {fj}, set z/ = in (j3.2ip . On the left-hand side we then have t{i^i), which 
is nonsingular since t^v) is entire, multiplied by g(fi) which is zero by (13.221) . Thus the 
left-hand side vanishes, and rearranging we have 

qjuj - 2ir]) ^ _ _2t^ a^ii^i,i]) , . 

or, using (|3.22p one more time, 

r ..n A sinh(2ir? + vi) _ a,^ a^(z^i, r?) i - i „ (3 24) 

^ \Uin^2ir,-n + l^^)~ b^{u.,r,) ' i3.24j 
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These are exactly the Bethe ansatz equations (|3.15p for the problem, with the v.i the 
roots. The expression for t[v) implied by (j3.2ip then matches the formula (|3.14p . as 
would have been found from a direct application of the Bethe ansatz. 

3.6 The quantum Wronskian 

In this review, the TQ relation will be our main tool in making the link with the theory 
of ordinary differential equations. However, there are other sets of functional equations 
associated with the six-vertex model which can be equally important. In principle these 
can all be obtained first as operator equations (see for example [26,31,33,132]), and then 
turned into functional relations by specialising to individual eigenvectors. However a full 
discussion of this would take us too far afield, so instead we shall profit from the fact 
that we have already obtained the TQ relation from the lattice model, and give some 
indications as to why these further properties should hold. 

We continue to consider the model with general (non-zero) twist, and discuss an 
important consequence of the identity ()A.48P : 

to(i^, 4>) = to{i^, -<P) = to{i^, 101) , (3.25) 

where to is the ground state eigenvalue of T(i/, (p) (with N/2 even) and (j) is the twist 
parameter. In view of ()3.25p . the following two TQ relations hold simultaneously: 

toii^,\cl)\)qoii^,<P) = a^ii^,rj)e~'^qoii^ + 2ir],<P)+b^iiy,7])e'^qoii^-2ii],<P); (3.26) 
to{u,\(t>\)qo{iy,-4>) = a^(z.,r/)e'%(z. + 2i?7,-,^) +6V,^)e-'%(z^-2ir/,-</.)(3.27) 

where go is the corresponding ground state eigenvalue of the corresponding matrix Q(z^, (j)). 
As equations for q, these only differ in the way that the twist factors e^**^ appear, and 
even this difference can be eliminated by defining 

qo{^,cP):=e-'''^/('^\o{i^,^) (3.28) 

so thalEI 

to{iy,\(l>\)qo{'^,(t>) = a^{u,rj)qo{u + 2ir],^)+b^{u,ri)qo{u-2ir],^); (3.29) 
to{u, \<P\) qo{u, -<P) = a^{i^, 7])qo{^ + 2irj, -0) + 5^(i/, rj)qo{iy - 2ir/, -<P) . (3.30) 

Thus qo{i^, (p) &iid qo{i^, —0) both solve the single functional equation 

to{i^A(l>\)qiiy)=a^{iy,7])qiu + 2i7])+b''iu,7])q{iy-2ir]). (3.31) 

This is a finite-difference analogue of a second-order ordinary differential equation, and 
so it should have two linearly independent solutions; equations (j3.29p and (j3.30p confirm 
that this is indeed the case. The quasi-periodicity in v induced by the definition ()3.28p . 
combined with the periodicity of the 'potential' to, means that qQ{i^,(p) and qo{i^,—(j)) 
can be interpreted as the two Bloch-wave solutions to (13.3ip [31, 133]. Just as in the 
continuum case, given two solutions to a single second-order equation it is natural to 

''in the field theory context, go will correspond to the vacuum eigenvalue of the operator Bazhanov, 
Lukyanov and Zamolodchikov denote Q, while go is the vacuum eigenvalue of their A. 
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construct their Wronskian. To this end, we can multiply (|3.29|) by qo{i', —(/)) and (|3.30p 
by qoii', (p), subtract and regroup terms by defining 



A{u) := qo{u+ir], -<j))qQ{v-i'q, 0) - qo{u+ir], 4>)qQ{v-i7^, -(/)) (3.32) 

to find 

= a^{u, ri)A{u+iri) - b^{u, r])A{u-ir]) . (3.33) 

Recalling from (j3.5p the definitions a{u, rj) = sin(r/ + ii/), b{i', rf) = sm{r] — iv), and the fact 
that is even, (13.330 implies that the function yV{v) := A(i/)/ sinh^(z^) is periodic with 
period P = 2ir]. However, from (j3.22p and (j3.28p . we see that W(z^) also has the period 
P' = 2iii. For P' / P = rj/TT irrational, W{v) must therefore be constant; by continuity in 
r/, W{v) is constant for all values of ry. 

Evaluating W at ^ oo gives an identity, the finite-lattice version of the quantum 
Wronskian relation discussed by Bazhanov, Lukyanov and Zamolodchikov in [31]. Re- 
expressed in terms of qQ{y., (p) it reads 

e~'^^qoii'+i'n, (j))qo{i'-irj, -(p) - e^'^qo{i'+ir], -(p)qo{iy-iri, (p) = -2ism{(p) sinh^(z^). 

(3.34) 

In deriving this result we have only treated the ground state eigenvalues of T and Q, 
and we have also assumed that the twist cp is non-zero. If these restrictions are dropped, 
a number of important subtleties arise, particularly in the so-called 'root of unity' cases 
when P' / P is rational. For more extensive discussions which address some of these issues, 
see, for example, [134] and [135,136]. 



3.7 The fusion hierarchy and its truncation 



A further functional equation results if we multiply the simultaneous TQ relations (j3.29p 
and (|3.30p by qQ{i>—2ir], —cp) and qo{i'—2ir], cp), respectively, and subtract. We then obtain 
not the quantum Wronskian, but instead an alternative expression for to in terms of % : 



to{i^,\<P\) 



-1 



2i sin ( 



qQ{i^+2ir], -(p)qo{u-2ir], cp) - qo{i^+2ir], (p)qo{u-2iri, -<p 
qo{u, (p)qo{iy-2ir], -cp) - qo{i^, -(p)qo{iy-2ir], (p) 

qo{v+2ir],-(p)qQ{v-2ir],(j)) - qo{i^+2iri, (p)qo{i^-2iri, -<p)){3.35) 



using the previously-obtained formula for the quantum Wronskian for the last equality. 
This is a sign that the quantum Wronskian fits into a hierarchy of relations which we 
now describe. It is convenient to change the normalisations slightly. Defining 



1 



y/—2i s'm< 
with rj = —r]+7r/2 and 



-"'^/^qo{u-ik^,(P),e"'^/^qo{u-ik?j, 



W[k,-k]{i^) :=det(g(^),Q 



(g(^))i {q^-% 



we set 



.-}^[k + l^-k-l]{u) , A: = -1,0,1,2,3... 



(3.36) 

(3.37) 
(3.38) 
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Then t(-i/2)(i.) = 0, and (ICTD and (IH:^^ imply 

(u) = [i cosh(i/)]^ , ^ ^ (3 39) 

We can then use the following Pliicker-type relation 

det(oo, Si) det(6o, &i) = det(6o, oi) det(ao) ^i) + det(6i, oi) det(6o) ao) (3.40) 
and the property 

W[k + a,-k + a]{ij) = W[k, -k]{u - ia?j) (3.41) 

to show that 

(i/ - i(2m+l)^)t(°) (u + i(2m+l)^) + ^("-1/2) (^)^(m+i/2) ^ (3^42) 

where the index m takes the half-integer values 1/2, 1, 3/2, .... Another set of relations 
among the functions t^'^^v) is also a simple consequence of the identity p.40p : 

i(i/2)(j,)iM(^_^(2m+l)^ = 

t(0)(j,_i^)i("^+i/2)(j,_^2m^) +t(°)(zy + i^t(™-i/2)(^_ •(2^^2)^) . (3.43) 

The sets of functional relations (|3.42p and (|3.43p are called fusion hierarchies [31,137-140]. 
The name comes from the fact that they can also be obtained by a process known 
as 'fusion' of the basic transfer matrix T, without introducing the auxiliary function 
q{v) [141]. 

An important phenomenon occurs at rational values of ry/vr, known as truncation of 
the fusion hierarchy. Here we shall just mention the case rf = 2^+2 ^ ~ 2M+2 ) '^'^^^ 
2M e Z+ and (p = 2m+2 - -^^^ periodicity of q{y) we have 

t(A^+V2)(^) = (3.44) 

and also 

( in iir 

which on comparing with 



(3.45) 



^ V-2isin</. \ ''0'^ 2M+2''>^^''^ WV*^ 2Af+2 ' 

( in in 
e4^+4g0(^+2fe>'A)'e"^^+^90(^^+2fe, 



(3.46) 



shows that 

i(*^)(z.) = det(g(2A/+i)^g-.(-2M-i)) ^ det(g(i),g(-i)) = tW(z.) . (3.47) 

Thus the infinite fusion hierarchy has been reduced, or truncated, to a finite set of func- 
tional equations (a T-system) constraining the t-functions. It is also easy to check that 
the relations (j3.42p . (j3.44p and (13.471) together imply the symmetry 

(u) = (u) , m = 0, 1/2, . . . M/2 . (3.48) 
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Truncation is important because it leads to closed sets of functional relations. Subject 
to suitable analyticity properties, these can be converted into sets of integral equations 
which allow the model to be solved - an example of this procedure will be given in 
appendix lD.lt in the slightly-simpler context of the large- limit. At generic values of (j) 
and rational ry/vr, truncation is also possible, though it has to be implemented in a more 
complicated way to ensure that the analyticity properties required for the derivation of 
the integral equations continue to hold [31,141,142]. 



3.8 Continuum limit of lattice models 

Very often, physicists are particularly interested in the behaviour of lattice models in the 
so-called 'thermodynamic limit', when the size of the system tends to infinity. If this 
limit is taken in a suitable way near to a phase transition, short-distance details of the 
model get washed away. The resulting behaviour is then said to be 'universal', and since 
it does not depend on the precise formulation of the model, it also tells us about more 
realistic systems near to their phase transitions, beyond the idealised models discussed 
so far. If we continue to measure distances by the number of lattice sites, all of the 
universal properties will be found in the long-distance asymptotics of quantities such as 
the transfer matrix eigenvalues to- To focus on these features, it is common to introduce 
a dimensionful lattice spacing d - up to now this has been equal to one - and then let this 
spacing tend to zero while keeping the 'physical' width of the lattice, Nd, finite. This 
process - which may have to be accompanied by a suitable tuning, or renormalisation, of 
parameters to ensure that the objects of interest retain finite values - is known as taking 
the continuum limit. It has the additional feature that the limiting theory can often be 
studied using techniques from quantum field theory. 

The six-vertex models lie at a phase transition of the more general eight- vertex model 
for all < T] < 7r/2, and so are well-suited to the taking of this limit. One place where 
universal behaviour can then be detected is in the behaviour of the logarithm of the 
dominant eigenvalue of the transfer matrix, to. As — > oo , 

lntoiN) = -fN+^ + ... . (3.49) 

The constant / is simply the large- limit of the free energy per site (|3.9p . and such a 
term is expected on general grounds. The next term is the signal of a phase transition: it 
depends only algebraically on the system size, and its general form is a consequence of the 
scaling symmetry characteristic of (second-order) phase transitions. If we now introduce 
the lattice spacing d, replace by L := Nd and define the subtracted/rescaled free 
energy to be 

F := -lnto{L) - fL (3.50) 

then the '. . .' terms in (I3.49P give vanishing contributions as c? ^ with L held fixed and 
in this limit 

nL) = -^, (3.51) 

where L is now a continuous (positive) number. This is the expected behaviour of the 
free energy for a conformal field theory (CFT) on an infinite cylinder with circumference 
L. In unitary theories with periodic boundary conditions, the proportionality constant 
Ceff coincides with the standard Virasoro conformal central charge c. The continuum 
limit of the six-vertex model and the XX Z spin chain is described by a unitary CFT 
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with central charge c^^ = 1. The effective central charge is c^g = c^^ = 1 in the periodic 
case, or 

'=."-^-^(^<' 

in the twisted cases [143,144]. 

The eigenvector ^^^^ corresponding to the dominant eigenvalue to just discussed is 
called the ground state of the model, and most of the results to be described later concern 
the Bethe roots for this state. However it is worth noting that the conformal field theory 
dictates the behaviour of all states in the continuum limit of the model [145,146]. The 
remaining states are sometimes called 'excited states'; they can be assigned an 'energy' 
as minus the rescaled logarithm of the corresponding eigenvalue, just as was done for 
the ground state to obtain the free energy F. The spectrum of the six-vertex model or 
XX Z model with periodic boundary conditions then consists of states with energies that 
behave as [128, 146-149] 

^{m,},{m[}Ak')iL) = ^ (^-^^ + ^i^k,k' +rn + m')^ , (3.53) 

where k, k' € Z, m = ^irrii, m! = '^^m'^ with rrii and m'^ non-negative integers, and 
Xk,k' = k'^^ + fc'^/(4x) with X = {tt — 2r])/2TT. The quantity ^ is a model-dependent 
parameter (the 'velocity of light'). It is 1 for the six- vertex model and, in the notation 
used here, tt sm2r]/2r] for the XX Z model. (Alternatively, it could be made equal to 1 
simply by multiplying the XX Z Hamiltonian by an overall factor [150].) For each pair 
k, k' , the state with m = m' = is associated in the CFT with what is called a 'primary 
field' with scaling dimension x. The states with other values of m and m' are called 
'descendants' of this field, and the even spacing of the energies of these descendants is a 
characteristic feature of the spectrum of a CFT. A similar tower-like structure also arises 
when twisted boundary conditions are imposed [147]. The terminology of primary fields 
and descendants relates to an underlying symmetry of the conformal field theories, the 
infinite-dimensional Virasoro algebra. For more on these topics, see, for example, [34]. 

It turns out that the universal behaviour described above corresponds to a special 
limit of the TQ and Bethe ansatz equations, in which their forms simplify. For this, it 
will be more convenient to re-express the results obtained in sections 13.21 and 13.51 using 
an alternative set of variables. Setting 

£^ = e2^' , a; = -6-2^" = e^^^ , (3.54) 

the BAE become 



n 

1=1 



N 



E'l - UJ^E'- \ _ 2n-N -i2^ ( ^+^E[ i - 1 n fS 55) 



If we continue to concentrate on the ground state, then, as already mentioned, n = A^/2, 
and all of the Vi lie on the real axis. This translates into all of the E\ being real and 
positive, and eliminates the factor w^"""^ from the RHS of the BAE. In the limit — > oo, 
the number of roots needed to describe the ground state diverges. This complication is to 
some extent compensated by the fact that the Bethe ansatz equations for the 'extremal' 
roots, those vi lying to the furthest left or right along the real axis, simplify in this limit. 
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at least for r] > tt/A. Since the left and right sets of extremal roots are constrained by 
the symmetry 

qo{-u, (j)) = qo{v, -(/>) ^ = -T^N/2+i~i{-<f) , (3-56) 

without loss of generality we shall concentrate on the left edge only. 

The left edge of the root distribution tends to — oo, as — ^logA^. (This behaviour 
can be extracted, for example, from the results of [28].) Hence the lowest-lying E'^ scale 
to zero as 

Ei ~ EiN-^i/"" . (3.57) 

To capture their behaviour, replace each E'- with N~^^/'^Ei in the BAE, and then hold 
the Ei finite as N ^ oo. The BAE simplify to the following: 



1=1 



(For r] < IT /A, the product must be regulated to ensure convergence in the N ^ oo limit, 
which complicates the story. We won't discuss this any further here, except to remark 
that it corresponds to the leaving of the semiclassical domain mentioned below.) 

A similar limit can be performed on qo^i'), taking care to adjust its normalisation to 
ensure a finite and non-zero result as — > oo : 

Qoiu) ^ qo{E) := lim le^^/^oii^)] , = ff f 1 - ^ V (3.59) 

Taking the same limit with tQ{u) and in general on the functions t^'^\v) results in a 
simplification of the TQ relation to 

to{E)qo{E) = e''^qo{uJ^E) + e-'^ qo{uj-^ E) , (3.60) 

and the fusion relations (j3.42p and (j3.43p become 

(L^"l^)t(") (loE) = 1 + V2) (^)^(m+l/2) (3_g^) 

and 

^(1/2) (^)iM (^2m+l^) ^ ^(^+1/2) (^2m_g^ ^ ^(m-1/2) (^2m+2^) _ (3_g2) 

These equations control the distribution of the extremal roots in the large- limit. 
Physically, they are important because they turn out to determine the constant Cgfr which 
controls the leading finite- A^ corrections to the ground state energy - see app endix lEl and 
for example, [28]. 

As already mentioned, for lo a root of unity the fusion relations truncate. For r] = 
7rM/(2M+2) with 2M E Z+, and c/) = 7r/(2M + 2), the truncated set of t-equations can 
elegantly be written as: 

ih-l)/2 ^ 

t^'") (L^"i^)t(") (loE) = 1 + Yl (i(j') (^)) , m = 1/2, 1, . . . , (/i - l)/2 

i=i/2 

(3.63) 

where h = 2M, uj = e'^'/(^^+^) and Gab is the incidence matrix of the Ay^^i Dynkin 
diagram: 
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1 



2 



3 



n-1 



n 



In the case of M = 3/2 and M = 2 the equations are respectively 



(3.64) 



and 



l + (t(l/2)(^))2 , 



(3.65) 



A first hint of a link with the theory of ordinary differential equations comes on comparing 
equation (j3.64p with the title and content of Sibuya's paper [151]: 'On the functional 
equation /(A) + f{u}X)f{u}~^X) = 1, (w^ = 1)'. A precise ODE/IM equivalence was first 
established in [1] by mapping (|3.65p into a functional relation which had previously been 
associated with the quartic anharmonic oscillator by Voros [21], and more importantly 
by showing an exact equivalence between the functions - and not just the functional 
relations - involved in the two setups. As explained in appendix [El techniques developed 
in the study of integrable models allow functional equations of the type (j3.63p to be 
transformed into sets of nonlinear integral equations known as thermodynamic Bethe 
ansatz (TBA) equations [152]. 

In conclusion: starting from the six-vertex model with twisted boundary conditions 
and using the algebraic Bethe ansatz approach we have derived sets of functional relations: 
the Baxter's TQ relation, the quantum Wronskian and the fusion hierarchy. Anticipating 
the correspondence with the theory of ordinary differential equations the continuum limit 
was taken, with the final set of equations encoding information about the c = 1 conformal 
field theory defined on a cylinder with twisted boundary conditions, with the value of 
the twist depending on the twist in the original six-vertex model. 

The precise way to extract information such as the effective central charge Ccg and the 
scaling dimensions goes through the transformation of functional relations into nonlinear 
integral equations (see appendix ID. ip . However, it turns out that there are other ways 
to derive the same sets of functional and integral equations rather than starting from 
the six-vertex model. One possibility is to work directly in field theory and exploit the 
fact that c = 1 CFT also corresponds to the ultraviolet limit of the sine-Gordon model. 
The derivation of the integral equations makes use of the sine-Gordon scattering matrix 
description and as mentioned before goes under the name of the thermodynamic Bethe 
ansatz [152]. 

Another approach also directly based on a CFT was proposed by Bazhanov, Lukyanov 
and Zamolodchikov in [30] and further developed in [31,32]. The starting point of [30-32] 
is not the unitary c = 1 conformal field theory defined on a strip geometry with different 
boundary conditions as above, but a CFT with central charge 



with periodic boundary conditions. This theory is neither unitary nor minimal and at 
fixed values of /3 the Hilbert space still depends on a free parameter p A brief summary 
of results relevant for the ODE/IM correspondence is reported in the next section. 



c = 1 - (/? - /3-i)2 < 1 



< /? < 1 



(3.66) 
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3.9 TQ equations in continuum CFT: the BLZ approach 

In [30-32], Bazhanov, Lukyanov and Zamolodchikov showed how for integrable models 
the structures such as Baxter's T and Q matrices may also be studied directly using 
field-theoretic methods. BLZ considered a CFT with central charge c parameterised in 
terms of /3 according to (j3.66p . The description of the conformal spectrum of the XX Z 
model in terms of towers of states, each tower consisting of a highest-weight state \p) 
and its descendents, applies to all CFTs. In this case, the highest-weight states have 
conformal dimension Ap = (p//?)^ + (c — l)/24, where p is a continuous parameter. 

For each tower of states, BLZ define a continuum analogue of the lattice transfer 
matrix T, an operator- valued entire function T(s,p). In analogy with (j3.28p . they also 
define a pair of operator- valued functions Q±(s,p). Together, these operators mutually 
commute and satisfy a TQ relation 

T(s,p)Q±(s,p) = Q±(g2s,p) + Q±(g-2s,p) (3.67) 

with q = exp(i7r/3^). Within each tower the highest-weight eigenvalues 

T{s,p) = (p|T(s,p)|p) (3.68) 

and 

Q±(s,p) = (p|5^p/^'q±(s,p)|p) , (3.69) 

satisfy the TQ relation 

T{s)Q±{s) = e^^^'PQ±{q-^s) + e^^^'PQ±{qh) , (3.70) 

where P is an operator such that P|p) = p\p)- If we set rj = -1(1 — /3^), this relation 
between the six-vertex anistropy rj and the coupling /? ensures that = exp(2z7r/3^) is 
equal to = exp(— 4«7y) and the BLZ and continuum six-vertex TQ relations match 
perfectly. Moreover, for in the so-called semiclassical domain < < 1/2, the 
eigenvalue Q+{s) can be written as a convergent product over its zeros {s^}: 

CO y ^ 

Q+(s) = n 1-- , (Q+(0) = 1). (3.71) 
fc=i ^ 

Thus a set of Bethe ansatz equations follows from the TQ relation and the entirety in s 
of the eigenvalues: 




(3.72) 



The other elements of the lattice picture also appear directly in the continuum con- 
text. From the identity operator To(s) and Ti(s) = T(s), an infinite set of mutually 
commuting operators are built using the fusion relations: 

T,{qs)T,{q-'s) = l + T^_^{s)T^^^{s) , j = i,l,|.... (3.73) 

At rational values of the parameter the hierarchy truncates to a finite set of operators, 
just as in the lattice case. Alternatively, the operators Tj(s) are given directly in terms 
of the Q's: 

2isin(27rP)Tj-(s) = Q+{q^^+^s)Q^{q-^^-h) - Q+{q-^^-h)Q^{q^^+h) . (3.74) 
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Evaluating on the state \p) with j = 0, we find the continuum version of the quantum 
Wronskian 

Q+{qs)Q_{q-h) - q~^ Q +{q~^ s)Q _{qs) = 2isin(27rf)) , (3.75) 

from which we deduce Q-{s,p) = Q+{s, ~p). 

Since much of the current interest in the functional-relations approach to integrable 
models is focused on the continuum field theory applications, we concentrate on the 
above version of the TQ relation in the remainder of these notes. However it should 
be remembered that the link with the theory of ordinary differential equations applies 
equally to lattice models, so long as a large-iV limit is taken in a suitable way. To be 
more precise, by setting 

^^ = l-^,P=f (3.76) 

one gets a full match between the six-vertex and the BLZ functional relations. Fur- 
thermore, for a non-unitary CFT on a cylinder with periodic boundary conditions the 
effective central charge is given by 

Ceff = C - 24Amin ■ (3.77) 

We therefore see that = c — 24Ap=o = 1 corresponds to the effective central charge 
for the untwisted six-vertex model. And, within a single p sector, the effective central 
charge associated to the highest-weight state \p) is 

c2=c-24Ap = l-24(^^j , (3.78) 

which matches the effective central charge for the six- vertex model with twisted boundary 
conditions (p. 

Bazhanov, Lukyanov and Zamolodchikov also discovered a relationship between the 
T and Q operators and (perturbed) boundary conformal field theory [30], which was 
followed up in later work [153-156]. In the long term this probably corresponds to the 
most fertile ground for applications of the ODE/IM correspondence. But as a description 
of this would take us too far from the themes of this review, we address the interested 
readers to the papers [30,154-156]. 



3.10 Summary 

We conclude our survey of integrable lattice models with a summary of the vocabulary 
introduced thus far, specialising to the simplified case of the continuum limit. In order 
to solve the six-vertex model, it suffices to diagonalise the 

transfer matrix T 
which depends on the 

spectral parameter u 
and has 

eigenvalues f(i/) . 
These can be given in terms of the 

Bethe roots {Ei} 
which solve 
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Bethe ansatz equations H {^ Ei-^-^e^ = -e"'^"^ 
and these can be neatly encapsulated in Baxter's 

TQ relation t{v)q{v) = e~''t'q{u + 2ir)) + e''t'q{iy - 2ir)) . 
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4 Ordinary differential equations and functional relations 



Surprisingly, the functional equations found in the last section also govern the problems 
in "PT-symmetric quantum mechanics discussed in section [21 To understand how this 
comes about, we must first return to the subject of "PT-symmetric eigenvalue problems 
and their generalisations in a little more depth. 



4.1 General eigenvalue problems in the complex plane 

We begin with one piece of unfinished business from section [2j what goes wrong with 
the Bender-Boettcher problem at M=2, and what can be done to resolve it? In figures 
[2] and [3l the energy levels continued smoothly past M=2, but in fact this can only 
be achieved by implementing a suitable distortion of the problem as originally posed. 
Consider the situation precisely at M=2 : the Hamiltonian is — x^, an 'upside-down' 
quartic oscillator, and a simple WKB analysis (about which more shortly) shows, instead 
of the exponential growth or decay more generally found, wavefunctions behaving as 

exp(ibix'^/3) as x tends to plus or minus infinity along the real axis. All solutions thus 
decay, albeit algebraically, and this complicates matters significantly. The problem moves 
from what is called the limit-point to the limit-circle case (see [45,46]), and additional 
boundary conditions should be imposed at infinity if the spectrum is to be discrete. 

While interesting in its own right, this is clearly not the right eigenproblem if we wish 
to find a smooth continuation from the region M < 2. Instead, it is necessary to enlarge 
the perspective and treat x as a genuinely complex variable. This has been discussed by 
many authors, and is particularly emphasised in the book by Sibuya [20] , though the 
treatment which follows is perhaps closer to that of [22,23]. 

The key is to examine the behaviour of solutions as |x| — > oo along a general ray 
in the complex plane, in spite of the fact that the only rays involved in the problem as 
initially posed were the positive and negative real axes. The WKB approximation tells 
us that 

i^ix) ~ P(x)"i/^ e± rVmdt (4,1) 

as |x| —> oo, with P{x) = —{ix)'^'^ +l{l+l)x~'^ — E . (This is easily derived by substituting 
ip{x) = f{x)e^^^^ into the ODE.) Since the problem was set up with a branch cut running 
up the positive-imaginary axis, it is natural to define general rays in the complex plane 
by setting x = pe^^/i with p real, as illustrated in figure [TTJ 

For M > 1, the leading asymptotic predicted by (|4.ip is not changed if P{x) is replaced 
by — (ix)^^^, and substituting into the WKB formula we see two possible behaviours, as 
expected of a second-order ODE: 

V'i ~ exp [±^e'^(^+^V'+''] • (4.2) 

For most values of 9, one of these solutions will be exponentially growing, the other 
exponentially decaying. But whenever 5Re[e*^(^+^^)] = 0, the two solutions swap roles and 
there is a moment when both oscillate, and neither dominates the other. The relevant 
values of 9 are 

" ■ " ■ ■ ... (4.3) 



2M+2 2M+2 2M+2 

(Confusingly, the rays that these values of 6 define are sometimes called 'anti-Stokes 
lines', and sometimes 'Stokes lines'. See, for example, [157].) 
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Figure 11: A ray in the complex x-plane. 



Whenever one of these lines lies along the positive or negative real axis, the eigenvalue 
problem as originally stated becomes much more delicate. Increasing M from 1, the first 
time that this happens is at M = 2, the case of the upside-down quartic potential 
discussed above. But now we see that the problem arose because the line along which 
the wavefunction was being considered, namely the real axis, happened to coincide with 
an anti-Stokes lin^. We also see how the problem can be averted. Since all functions 
involved are analytic, there is nothing to stop us from examining the wavefunction along 
some other contour in the complex plane. In particular, before M reaches 2, the two ends 
of the contour can be bent downwards from the real axis without changing the spectrum, 
so long as their asymptotic directions do not cross any anti-Stokes lines in the process. 
Having thus distorted the original problem, M can be increased through 2 without any 
difficulties. The situation for M just bigger than 2 is illustrated in figure \T2\ with the 
anti-Stokes lines shown dashed and the wiggly line a curve along which the wavefunction 
ip{x) can be defined. 
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Figure 12: A possible wavefunction contour for M > 2. 
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as just mentioned, some would have called this a Stokes line. 
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The wedges between the dashed hnes are called Stokes sectors, and in directions out 
to infinity which lie inside these sectors, wavefunctions either grow or decay exponen- 
tially, leading to eigenvalue problems with straightforward, and discrete, spectra. Note 
that once M has passed through 2, as in figure [T2l the real axis is once again a 'good' 
quantisation contour - but for a different eigenvalue problem, which is not the analytic 
continuation of the original M < 2 problem to that value of M. (For the analogue of 
figure E] for this new problem, see figure 20 of [24].) Going further, we could choose any 
pair of Stokes sectors for the start and finish of our contour. A priori, each pair of sectors 
defines a different problem, though we shall see later that some of these problems are 
related by simple variable changes. 

All of the problems do share one feature - their quantisation contours begin and end 
in the neighbourhood of the point x = oo. In the terminology of the WKB method, they 
are related to 'lateral' connection problems [158]. There is one other special point for 
the ordinary differential equation, namely the origin, and this provides another natural 
place where quantisation contours can end. Contours which join x = to x = oo lead to 
what are called 'radial' (or 'central') connection problems, and with suitable boundary 
conditions they can also have interesting, discrete, spectra. However, if both ends of the 
contour are placed at x = 0, the resulting eigenvalue problem is always trivial. Some 
sample quantisation contours are shown in figure [T3l 



We should pause for a moment to consider which boundary conditions can be imposed 
at the origin, in order to understand why x = and |x| = oo behave differently as end 
points of quantisation contours. Even with the angular-momentum-like term l{l+l)x~'^ 
included, the singularity at the origin is much milder than that at |x| = oo, and irre- 
spective of the direction in which it is approached, solutions there behave algebraically, 
as x'+^ or x~'. For this reason the complications associated with Stokes sectors do not 
arise in the neighbourhood of the origin, and there are just two natural boundary con- 
ditions to impose there - we can either demand that the solution behaves as x'+^, or 
as x~^ (the more singular of these two boundary conditions being defined by analytic 
continuation). This contrasts with the situation near |x| = oo, where we can ask that 
a solution be subdominant in any one of the potentially infinitely-many different Stokes 
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Figure 13: Some further quantisation contours. 
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sectors. Expressed more technically, the ordinary differential equation has two singular 
points, one at the origin and one at infinity. The singular point at the origin is regular, 
and solutions have straightforward series expansions in its vicinity. These converge in 
the full neighbourhood of the origin, and can be analytically continued in a simple 
Infinity, on the other hand, is an irregular singular point, in the neighbourhood of which 
solutions have asymptotic expansions which only hold in selected Stokes sectors. This 
makes analytic continuation around the point at infinity much more subtle, and indeed 
this will be a major theme in the subsequent development. 

To summarize: associated with an ODE of the type under consideration there are 
many natural eigenvalue problems, which fall into two classes. Problems in the first, 
lateral, class are defined by specifying a pair of Stokes sectors at infinity, and then 
asking for the values of E at which there exist solutions to the equation which decay 
exponentially in both sectors simultaneously. Problems in the second, radial, class are 
defined by demanding decay in a single Stokes sector at infinity, and imposing one of the 
two simple boundary conditions at the origin. The questions in PT-symmetric quantum 
mechanics discussed in section [2] are all related to lateral problems, with one particular 
pair of Stokes sectors selected. Considerations of analytic continuation have led us to 
put all pairs of sectors on an equal footing, and we completed the story by bringing in 
the radial problems as well. But at this stage each eigenvalue problem sits on an isolated 
island, each with its own private spectrum. 

In the next subsection, we shall start to construct some bridges between the islands, 
using methods inspired by earlier work of Sibuya and of Voros. Remarkably, these bridges 
turn out to be precisely the functional equations which had previously arisen in the 
context of integrable quantum field theory. 



4.2 A simple example 

To illustrate the basic ideas in the simplest possible way, for the time being we set 
/(/+!) = 0, so we are dealing with the original Bender-Boettcher family of eigenproblems: 

- T^'iPix) - {ixf^^i;{x) = Ei}{x) , V G L^{C) ■ (4.4) 

Now that the perspective has been widened to encompass eigenvalue problems on 
general contours, it is convenient to eliminate the factors of i appearing everywhere by 
making the variable changes 

x^x/i, E ^ -E (4.5) 
so that the differential equation becomes 

tp{x)={). (4.6) 

This also moves the branch cut onto the negative real axis, and the initial quantisation 
contour onto the imaginary axis. (This final point explains why the reality of the spec- 
trum remains a non-trivial question, despite the fact that the coefficients in (14. 6p are all 
real.) 

^Cases with 2M not an integer fall just outside the treatments in the standard texts, but they behave 
in essentially the same way - see, for example, [159]. We have also glossed over some details, such as the 
logarithms which can arise in certain situations. More background on these issues can be found in [160]. 



dx'^ 



+ x 



2M 



E 
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Next, we need to develop our treatment of ordinary differential equations in the 
complex domain, relying largely on the work of Sibuya and co-workers [20, 161]. The key 
result is the following: 

• The ODE (j4.6p has a 'basic' solution y{x, E) such that 

(i) y is an entire function of x and E; 

[Though, because of the multivalued potential, x lives on a cover of C\{0} if 2Af^Z0.] 

(ii) as oo with | arg x \ < 37r/(2M+2), 



/2i 
1 



X exp 



:x^^/2exp 



M+l 



Af+1 



(4.7) 
(4.8) 



[Though there are small modifications for Af < 1 - see, for example, [16] . ] 
Furthermore, properties (i) and (ii) fix y uniquely. 

The second property can be understood via the WKB discussion of section WTH With 
the shift from x to x/i, the anti-Stokes lines for (j4.6p are 



arg(x) 



vr 



Stt 



"2M+2 ' 2M+2 ' 
and in between them lie the Stokes sectors, which we label by defining 



arg(2;) 



27rk 



2M+2 



< 



TT 



2M+2 



(4.9) 



(4.10) 



Three of these sectors are shown in the figure a 90° rotation of figure [T2l 

The asymptotic given as property (ii) then matches the result of a WKB calculation 
in 5_i U5o U5i . The determination of the large behaviour of the particular solution 
y{x, E) beyond these three sectors is a non-trivial matter, since the continuation of a 
limit is not necessarily the same as the limit of a continuation. This subtlety is related to 
the so-called Stokes phenomenon, and it can be handled using objects known as Stokes 
multipliers, to be introduced shortly. 

One more piece of terminology: an exponentially-growing solution in a given sector 
is called dominant (in that sector); one which decays there is called subdominant. It is 
easy to check that y as defined above is subdominant in Sq, and dominant in 5_i and Si. 
A subdominant solution to a second-order ODE in a sector is unique up to a constant 
multiplier; this is why the quoted asymptotics are enough to pin down y uniquely. 

Having identified one solution to the ODE, we can now generate a whole family using 
a trick due to Sibuya. Consider the function y{x, E) := y{ax, E) for some (fixed) a G C. 
From (j4.6p . y{x,E) satisfies 



dx^ 



+ a^^+^x^^ 



a^E 



y{x,E)=0. 



(This is sometimes given the rather-grand name of 'Symanzik rescaling'.) 
shifting E to a~'^E shows that y{x,a~'^E) again solves ()4.6p . Defining 



UJ 



J2-Ki/{2M+2) 



(4.11) 
If a2A^+2=l, 



(4.12) 



^°The original work of Hseih and Sibuya [161] concerned only the case 2M £ N, but the result also 
holds for the more general situation 2M £ of eq. (|4.6|l , so long as the branching at the origin is taken 
into account. This generalisation was explicitly discussed by Tabara in [162]. 
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Figure 14: Three Stokes sectors for the ODE (jM]), with M = 2.1. 



and 

E) := u;'=/2y(a;-^x, w^'^i?) (4.13) 

we then have the statements 

• Vk solves (gS]) for all A; G Z; [This follows since {uj-^Y'^''+^ = 1.] 

• up to a constant, yk is the unique solution to ()4.6p subdominant in Sk- [This follows 
easily from the asymptotic of y. ] 

• the functions yk-, yk+i are linearly independent for all k, so each pair {yk.,yk+i} forms 
a basis of solutions for (j4.6p . [This follows on comparing the asymptotics of yk and yk+i in 
either Sk or Sk+i-] 

We have almost arrived at the TQ relation. Next, the fact that can be expanded 
in the {yQ-,yi} basis shows that a relation of the following form must hold: 



E) = C{E)yo{x, E) + C{E)yi{x, E) . 



(4.14) 



This is an example of a Stokes relation, with the coefficients C{E) and C{E) Stokes 
multipliers. They can be expressed in terms of Wronskians, where [163] the Wronskian 
of two functions / and g is 

W[f,g\:= fg' - fg. (4.15) 

Given two solutions / and 3 of a second-order ODE with vanishing first-derivative term, 
their Wronskian VF[/,(7] is independent of x, and vanishes if and only if / and g are 
proportional. As a convenient notation we set 



Wk^,k2 ■= W[yk„yk2] = VkiV'ki - y'kiVki 
and record the following two useful propertieJ^: 

Wk,+i,k,+i{E) = Wk,,k,{oo^E) , Wo,i{E) = 1 . 



(4.16) 



(4.17) 



Now 'taking Wronskians' of the Stokes relation (j4.14jl first with yi and then with yo 
shows that 

■ (4.18) 



c 



Wo. 



c 



Wo. 



1 



^^The normalisations in (|4.7|) and (I4.13|) . which differ from those adopted by Sibuya, were expressly 
chosen so as to simplify these two formulae. 



33 



and so the relation itself can be rewritten as 

C{E)yoix, E) = y^iix, E) + E) , (4.19) 
or, in terms of the original function y, as 

C{E)y{x, E) = io~^/^y{uJX, uj~^E) + Lo'^/^y{Lu'^x, lo^E) . (4.20) 

This looks very like a TQ relation. The only fly in the ointment is the a;-dependence of 
the function y. But this is easily fixed: we just set x to zero. We can also take a derivative 
with respect to x before setting x to zero, which swaps the phase factors a;±i/2. So we 
define 

D_{E):=y{0,E) , D+{E) := y' {0, E) . (4.21) 
Then the Stokes relation (I4.20p implies 

C{E)D^{E) = LO^^/^D^{uj-^E) + uj^^/'^D^ioj'^E) , (4.22) 

and precisely matches the forms of the TQ equations ()3.60p and ()3.70p if the twist pa- 
rameter is set io (f) = 2ttp = ■k/{2M+2). Although the details are at this stage sketchy 
- more will be provided in later sections - we can already see how some concepts in the 
two worlds of integrable models and ordinary differential equations must be related: 



Six-vertex model with twist 


Schrodinger equation with 


= 27rp = ■k/{2M+2) 


homogeneous potential x'^^ 


Spectral parameter 4- 


Energy 


Anisotropy 4- 


-> Degree of potential 


Transfer matrix 


-> The Stokes multiplier C 


Q operator 4- 


-> D^{E) : the value of y(x, i?) at x = 



If y on the last line is replaced by y\ then the twist changes io (f) = — 7r/(2M-|-2). 

A small puzzle remains at this stage: why should one particular value of the twist in 
the integrable model be singled out for a link with an ordinary differential equation? This 
was resolved shortly after the original observation in [1], when Bazhanov, Lukyanov and 
Zamolodchikov [2] pointed out that including an angular-momentum-like term l{l+l)/x'^ 
allowed Q operators at other values of the twist to be matched. The details, and a further 
small generalisation of the basic ODE (14.61) . will be covered in section [5l 

4.3 The spectral interpretation 

How should we think about C and Dl In fact they are spectral determinants. The 
spectral determinant of an eigenvalue problem is a function which vanishes exactly at 
the eigenvalues of that problem: it generalises to infinite dimensions the characteristic 
polynomial det(M — AI) of a finite-dimensional matrix. Recall that C{E) is equal to 
the Wronskian W-i^i{E). Thus C{E) vanishes if and only if W[y-i,yi\ = 0, in other 
words if and only if E is such that y_i and yi are linearly dependent. In turn, this is 
true if and only if the ODE (jS.ip has a solution decaying in the two sectors 5_i and 5i 
simultaneously, which is exactly the lateral eigenvalue problem discussed in section [2l 
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modulo the trivial redefinitions of x and E. This is enough to deduce that, up to a 
factor of an entire function with no zeros, C{E) is the spectral determinant for the 
Bender-Boettcher problenJ^. Even this ambiguity can be eliminated, via Hadamard's 
factorisation theorem, once the growth properties of the functions involved have been 
checked; see [4] for details. By its definition, the zeros of D- (E) are the values of E at 
which the function y{x), vanishing at x = oo, also vanishes at x = 0. Likewise, the zeros 
of D^(E) are the values of E at which y{x) has zero first derivative at x = 0. Thus 
Dzp{E) are also spectral determinants. Note that the vanishing of D_ or Z?+ corresponds 
to there existing normalisable wave functions for the equation on the full real axis, with 
potential jxp^''^, which are odd or even, respectively, as illustrated in figures [T5] and [T6| 
this explains the labelling convention adopted earlier. 




Figure 15: |x|-potential, even sector: the ground state wavefunction. 




Figure 16: |x|-potential, odd sector: the first excited state wavefunction. 

This insight allows a gap in the correspondence to be filled. We mentioned in sec- 
tion [3?2] that while the TQ relation is very restrictive, it does not have a unique solution. 
So to claim that D-{E) is 'equal' to Aj^[\,p) begs the question: which ^+(A,p)? To 
answer, we first note that the radial problem for which D-{E) is the spectral determi- 
nant, in contrast to the lateral Bender-Boettcher problem, is self-adjoint, and so all of 

Since we performed a variable change in this section compared with the discussion in section [21 it is 
in fact C'{—E) which provides the spectral determinant for the Bender-Boettcher problem as originally 
formulated. 
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its eigenvalues are guaranteed to be real. Back in the integrable model, we mentioned 
previously that there is only one solution to the BAE with all roots real, so the question 
is answered: the relevant ^+(A) is that corresponding to the ground state in the spin-zero 
sector of the model. 
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5 Completing the dictionary 



5.1 Adding angular momentum 

We now restore the angular momentum term, and consider the differential equation 



dx"^ x"^ 



^{x) = E^{x). (5.1) 



This is the '7r/2-rotated' version of the generalised Bender-Boettcher problem (|2.5p . In 
the process of discussing this ODE we shall fill in some of the technical details skipped 
previously. The I = case was the subject of the first ODE/IM correspondence in [1], 
while the generalisation to / 7^ was introduced in [2]. The initial treatment below 
follows [4]. 

At the origin, solutions to (15. ip behave as a linear combination of rc'^^ and x^K As 
mentioned in section 14.11 a natural eigenproblem for this equation asks for values of E for 
which there is a solution that vanishes as x — > +00, and behaves as x^^^ as x — > 0. In the 
Stokes/WKB language this is a radial problem. For ?R.el > —1/2, the condition at the 
origin is equivalent to the demand that the usually-dominant behaviour there should 
be absent; outside this region, more care is needed, but the problem can be defined by 
analytic continuation. This issue will be discussed in more detail in section [53] below. 

As for the simple example above we first employ Sibuya's trick. From the uniquely- 
determined solution y{x,E,l), with large, positive x asymptotic [20] El 



y{x,E,l) ^ —^exp\^-^j , (5.2) 

we generate a set of functions 

yk{x,E,l) =Lu''/^y{io~''x, LO^'^E, I) , u; = , keZ (5.3) 

all of which solve (j5.ip . As before, any pair {yj., yk+i} forms a basis of the two-dimensional 
space of solutions, and so y_i can be written as a linear combination of yg and yi. 
Rearranging the expansion and using the properties (j4.16p . which continue to hold with 
the addition of the angular-momentum term, 

C{E, I) yo{x, E, I) = y_i(x, E, I) + y i(x, E, I) , (5.4) 

where Stokes multiplier C{E, I) again takes a simple form in the normalisations we have 
chosen: 

C{E,l) = W[y^,,yi]/W[yo,yi] = W[y.,,yi] ■ (5.5) 

There is one last complication: the angular-momentum term in (j5.ip means x cannot 
simply be set to zero in ()5.4p to find a TQ relation. Instead y should be projected onto 
another solution, defined via its asymptotics as x — > 0. Given that solutions near x = 



The result concerning the entirety of y proved by Sibuya (see Ref. [20]) at I — 0, 2M £ N also holds 
for the more general situation of eq. (jS-lf) and eq. (jG-lfl below (so long as the branching at the origin is 
again taken into account). The I — 0, 2M G case was discussed in [162], while the generalisation to a 
potential P{x)/x'^ with P{x) a polynomial in x was studied by MuUin [164], and more recently in [165]. 
It is also worth noting that with a change of variable it is possible to map (|5.1|l and (16.11) with / £ R, 
2M G onto particular cases of those treated in [164]. 
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behave as linear combinations of and x a solution 'il){x,E,l) can be defined by 
the requirement 

V'(x,S,0 ~2;'+^ + 0(x^+^). (5.6) 

This defines ^ uniquely provided Ke/ > —3/2. A second solution can be obtained from 
the first by noting that the differential equation - though not the boundary condition - 
is invariant under the analytic continuation I —l—l. As a result, ip{x, E, —l—l) also 
solves (15. ip . Near the origin, it behaves as , and so at generic values of / the pair of 
solutions 

i;+{x,E) := i^{x,E,l), 

il^-{x,E) := Tpix.E.-l-l) (5.7) 

are linearly independent. There are some subtleties to this procedure at isolated values 
of /, to which we shall return in section below. However, they do not affect the initial 
argument. In discussions of the radial Schrodinger equation (see, for example, Ch. 4 
of [166]), ijjj^ is sometimes called the regular solution, if 5Re/ > —1/2. 

We now take the Wronskian of both sides of (j5.4p with ^/>(x,i?, /) to find an x- 
independent equation: 

C{E, l)W[yo, ij]{E, I) = VF[y_i, V](i^, /) + W[yi,i^]{E, I) . (5.8) 

To relate the objects on the right-hand side of this equation back to VF[yO)V']) we first 
define a set of 'rotated' solutions by analogy with (j5.3p : 

V'fc(x, E, I) = u^/^^Piuj-^'x, io^'^E, I), kGZ. (5.9) 

These also solve (j5.ip , and a consideration of their behaviour as x — > shows that 

V'fc(x, E, I) = w~('+i/2)V(^, E, I) . (5.10) 

In addition, 

W[yk,A]{E,l) = u;'W[yiu;-'x,u;^''EJ),i;{iv-'x,u;^'E,l)] 

= W[y,i,]{uj^'E,l). (5.11) 

Combining these results, 

W[yk,^l^]iE,l)=J'+'/^^'W[y,i^]{u;^''E,l) (5.12) 

and so, setting 

D{E,l):=W[y,i;]{E,l), (5.13) 
the projected Stokes relation (j5.8p is 

C{E, l)D{E, I) = c<j~('+i/2)^(^-2^^ ^ J^+^/^^D{uj^E, I) . (5.14) 
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5.2 Matching TQ and CD relations 

Finally we are ready to make the precise connection with the TQ relation (j3.70p . As a 
shorthand, define 

D^{E):=W[y,^±]{E,l) (5.15) 

(so D^{E) = D{E,l) and D+{E) = D{E,-l-l)). Then (|5l^ taken at I and -l-l 
b GConiGS 

C{E, l)D^{E) = LU^^^'+^y^D^iio-^E) + a;±(2'+i)/2^^(^2^) _ ^^^-^g) 

If we set 

then the match between the general TQ relation (jS.TOp and the Stokes relation (j5.16p 
is perfect, with the following correspondences between objects from the IM and ODE 
worlds: 

T ^ C 

The mapping could also have been made onto the limiting form of the Bethe ansatz 
equations for the six-vertex model with twisted boundary conditions that was obtained 
in section 13. 5i In this case we have 

t ^ C (5.18) 

q ^ D , 

and, as mentioned before, the relationships between the anisotropy parameter ry and the 
twist (j) of the lattice model, and the parameters M and I appearing in the potential of 
the Schrodinger equation are 

vr M (2/ + l)7r 

The exact mapping between the functions appearing in these equations can be found 
by examining the functions &i E = Q and their asymptotic behaviour at E = co. Since 
the behaviour of D^{E) can be deduced from that of D_{E) = D{E,l), we only need 
the following results [4], which hold for M > 1: 

1. C and D are entire functions of E; 

2. The zeros of D are all real, and if / > —1/2 then they are all positive; 

3. The zeros of C are all real, and if —1 — M/2 < I < M/2 then they are all negative; 

4. If Af > 1 then the large-S asymptotic of D is 

In D{E J) ^i-^)^ ' 1^1 ^oo , |arg(-^)| < vr , (5.20) 

where /x = (M+1)/2M, and = - r(-/i)r(^ + i)//?? ; 

5. At £^ = 

1 r,l , 1 21 + 1 I 1 

DiE, Ol^^o = ^r(l + ^)(2M+2)5MT2+5 . (5.21) 
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6. The large- ii^ asymptotic implies that D(E, I) has orde^o equal to /i, which is strictly 
less than one for M > 1. Thus, Hadamard's factorisation applies in its simplest 
form and D{E, I) can be represented as: 

D{E,l) = D{0,l)f[(l-^) (5.22) 



k=l 



Property (i) follows from the definition of -D as a Wronskian, since all functions involved 
are entire functions of E. Properties (ii) and (iii) will be proven in section [6?2l while the 
proof of properties (iv) and (v) can be found in appendix [Bl 

The relevant analytical properties of T(s) and Q+{s) are given in [30,31]. For in 
the semiclassical domain: 

1. T{s,p) and Q±{s,p) are entire functions of s with an essential singularity at infinity 
on the real axis; 

2. All zeros of Q+{s,p) are real, and if 2p > — they are all strictly positive; 

3. All zeros of T{s,p) are real, and if |p| < 1/4 they are all negative; 

4. The large-s asymptotics are 



lnQ±(s,p) ~ ao(M+l)r(2^) 



2fi 



where // = 1/(2-2/32) and ao is as defined above; 
5. If s = 



oo , arg(— s) < vr , 

(5.23) 



g+(0) = 1 . (5.24) 

6. The large-s asymptotic implies that Q±{s) has order equal to 1/(2 — 2/3^), which 
is again strictly less than one for 0^ in the semiclassical domain and we can write 

oo ^ \ 

Q±(s)=n 1-f (5-25) 

fc=l ^ 

The restriction of /^^ to the domain < /J^ < 1/2 translates to the constraints 
7r/4 < ry < 7r/2 on the lattice model parameter ry, and we see that the point at 
which the factorised products have to be regularised coincide in the two cases. 
For convenience we have only considered 0^ inside the semiclassical domain; see 
[4,31,32] for discussions on the interesting case of /J^ > 1/2. 

Given (j5.17p . properties (i-iii) and (vi) match the equivalent statements for T{s,p) 
and Q+is,p) with E replaced by s. Noting that ^ = (M+1)/2M = 1/(2-2/?^), the 
asymptotic (iv) and the normalisations of D- and can be made to agree by setting 

s = vE , i; = (2M+2)-i/^r(^)-2 , = D^{0,2p/ (3^ -^y' . (5.26) 



^■'Technically, the order of an entire function f{z) is defined to be equal to the lower bound of all 
positive numbers B such that \f{z)\ — ©(e'^' ) as \z\ — > oo. 
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The precise result is [4] 



Q±is,P)^, = 7T^T(f,f-^),,,^,_,_, ; (5.27) 



■s 2p _ 1 - 



A/=/3-2-l 



5.3 The role of the fusion hierarchy 



ns,p)^, = C^(f,f-^)__ (5.28) 



M=/3-2-l 



In section [3771 another set of functional relations found in integrable models was described: 
the fusion hierarchy. Now that the TQ relation has been mapped onto a Stokes relation, 
it is natural to ask whether an analogue of the fusion hierarchy can also be found in the 
differential equation world, and it turns out that this is indeed possible [4]. 

Previously we examined the expansion of y_i in the basis {yo, yi}, but one can equally 
ask about the expansion of in any other basis, such as {yk+r-ijUk+r}'- 

yk^i = C^^'\k+r-i + C^k^Vk+r (5.29) 
A change of basis from {yk+r-i,yk+r} to {yk^i.,yk} can then be encoded in a 2 x 2 

(r) 

matrix as 



This matrix depends on E and /, but not x. The following properties are immediate: 

cM{£;,i) = C<:^i(w^£,i), (5.31) 



i ? . C«=f^;' -}]. ,5.32) 



(1) 

Further relations reflect the fact that the change from the basis {yfc+r+n-i, yfc+r+n} to 
{yfc+r-i,yfc+r}, followed by the change from {yk+r^i^Vk+r} to {yk-i,yk}, has the same 
effect as accomplishing the overall change in one go: 

cV^tl = Cr™^ . (5.33) 

(These express the consistency of the analytic continuations, and can be thought of as 
monodromy relations.) The case r = 1 gives two non-trivial relations: 

C^k^C^:^.-Ct2^=ct^'^ (5.34) 

and 

•^fc '^k+1 ~ ^k+2 - ^k ' 

which can be combined with the 'initial conditions' (|5.32p : to deduce cj^^ = —Cj^^ ; and 
then the more general equality 

) = -Cf-'^ (5.36) 
follows on comparing ()5.35p with ()5.34p . If we now set 

C7(") (E) = C7^") (a;-"+i£;) , (5.37) 
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then (|5.34|) is equivalent to 

C7(^)C(")(tj"+^^) = C(""^)(w"+2^) + C("+^)(l<j"^) , (5.38) 

and this matches the fusion relation (j3.62p . Since C(°)(S) = 1 = To{E) and, from the 
last section, C''^\E) = C{E) = Ti/2{vE), this establishes the basic equality 

C^''\e)=T^;,{vE) . (5.39) 

To find the fusion relation (|3.6ip . we start by taking Wronskians of (j5.29p as before 
to obtain 

C^"^ = Wk-l,k+r , C^k^ = -Wk-l,k+r-l , (5.40) 

from which we immediately recover (j5.36p . and can also deduce 

This relation combined with the n = —r case of (jOgi) . namely C^ 'C^^; = 1, implies 
that 

C(^-i) (cj-i^)C('-^) (ioE) - (S)C7(^-2) {E) = l , (5.42) 

which reproduces the fusion relations ()3.6ip . given the identification (j5.39p . 

Finally, combining (j5.39p and (|5.40p we have an expression for Tn/2ivE) in terms of 
a Wronskian: 

T,/2{vE) = C^-\E) = W^i,n{^-''^'E) . (5.43) 

This result shows that the fused transfer matrices can also be interpreted as spectral 
determinants. The right hand side of (j5.43p vanishes if and only if E is such that y_i and 
Un are linearly dependent, which in turn is true if and only if the ODE has a nontrivial 
solution which simultaneously decays to zero as \x\ — > oo in the sectors 5_i and S^. This 
is one of the lateral eigenvalue problems discussed in section 14.11 with the eigenvalues 
encoded in the zeros of C^"^ The full story is illustrated in figure [T71 
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Figure 17: General quantisation contours involved in the ODE/IM correspondence. 

Once the fused transfer matrices have been understood in this way, truncation of the 
fusion hierarchy can be reinterpreted in terms of the (quasi-) periodicity (in k) that the 
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functions exhibit whenever M is rational [4]. In the simplest cases (with M rational 
and /(/+1)=0) this periodicity arises because the solutions to the ODE live on a finite 
cover of C\{0}; for other cases, the monodromy around x={) needs a little more care, but 
the story remains essentially the same. 

As a simple example, consider the Schrodinger problems with 2M integer and /(/ + 
1) = 0. In these cases all solutions of the ODE are single- valued functions of x, and 
the sectors Sn+2M+2 and Sn coincide. Thus both y„ and yn+2M+2 are subdominant in 
Sn and must be proportional. In fact, one can easily show from the asymptotics that 
yn+2M+2{x^ E, I) = —yn{x, E, I), from which we conclude from (I5.43P that 

Cm)(^E) = l , C(2*^+^)(i?) = . (5.44) 

Thus the set (I5.42p of functional relations truncates to 

2M-1 

C^'\uj-^E)C^'-^ujE) = 1+ n {C^''HE)f"'^ , (5.45) 

n=l 

perfectly matching the T-system (|3.63p . 
5.4 The quantum Wronskians 

There is one final set of functional relations to discuss: the quantum Wronskian (j3.75p 
and its partner relations (j3.74p which express the T's in terms of the Q's. Recall the 
solutions ^/J± , defined in equation ()5.7p by their behaviour at x = 0. At generic values of 
I, these give an alternative basis in which to expand y{x, E, I). Using the Wronskian 

W[tl;+,^^]=2l + l (5.46) 

(evaluated in the limit x 0), and the relations Dzp = W[y,^p±] , 

{2l+l)y{x,E,l) = D^{E)^^{x,E) - D+{E)^+{x,E) . (5.47) 

Notice that there is a problem with this expansion at I = — 1/2 , which is easily under- 
stood: at this point the two solutions ■0+ and ■0- coincide and no longer provide a basis, 
a fact which is reflected in the vanishing of their Wronskian (j5.46p . In fact this is not the 
only place where difficulties arise. So long as I lies in the right half-plane 'Re I > —1/2, the 
solution tp+i^i = ip{x, E, I) can be proved to exist as the limit of a convergent sequence 
of approximate solutions - see [159,166,167]. However, this does not work in the left half- 
plane, which is another way to see why the second solution tp-{x,E) = il^ix^E,—!—!) 
must initially be defined by analytic continuation. At isolated values of I in the left half- 
plane, poles may arise, causing to be ill-defined. As discussed in Chapter 4 of [166], 
these poles can be removed by multiplying ■0_ by a regularising factor. However this 
inevitably inserts extra zeros into the Wronskian (j5.46p . and the regularised ip- fails to 
be independent of ■0+ at exactly the points where the previous ip- had failed to exist. 
For the simple power-law potential x^*^ that we are dealing with, the problem values of 
/ are best identified via the iterative construction of [159], and lead to the conclusion - 
see for example [4] - that {■0+,'!/'-} fails to be a basis at the points 

i + i = ±(mi + (Af+l)m2), (5.48) 
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where mi and m2 are two non-negative integers. Using (|5.17|) . for the integrable model 
this corresponds to the twist values 

2p = ±{mil3'^ + 1712} . (5.49) 

The set 

2p = ±m2, (m2 = 0,l,...) (5.50) 

corresponds to the vanishing-points for the quantum Wronskian (see eq. (I3.75P ) and at 

2p = -(3^-m2, (7712 = 0,1,...) (5.51) 

there is a normalisation problem for Q^{s^p) as a consequence of the appearance of a 
zero level {sk = for some k) [31]. 

The problem points (|5.48p can be dealt with by a limiting procedure - [4] discusses 
the first case, / = —1/2. For now, though, we will assume that I has been picked so that 
these subtleties do not arise. Then the pair {ipj^^ij:^} does provide a basis for the ODE, 
and, using (|4.13p . we define pairs of solutions V'fc^(x,£') as 

i;^{x,E) =uj''/^^±{uj-''x,oj^''E,l) , keZ. (5.52) 

These allow the rotated functions to be expanded as in ()5.47p : 

{2l+l)yk{x,E,l)=D^iuj'^''E)ij-{x,E) + D+{uj'^''E)i^^{x,E). (5.53) 

The Wronskians of these solutions are very simple: 

W[^^^^,^l:;] = W[ij^,,jj;] = 0, (5.54) 

= (2/+l)c^('^-P)('+i/2). (5.55) 

The fundamental quantum Wronskian is now almost immediate: in our normalisations, 
yo] = 1 ! and substituting the expansion (j5.53p into this formula and using bilin- 
earity of the Wronskian yields 

{21 + 1) = io~'^^+^/^^D^{u~^E)D+{ijE) - J+^/^D^{iuE)D+{io-^E) . (5.56) 

The fused Wronskians are equally straightforward. Taking the Wronskian of y_i and y„, 
again using the expansion (j5.53p . shifting E to w^^^^E and using the relation (j5.43p for 
C7(") shows that 

(2/ + l)C(")(^) = w-("+i)('+i/2)Z)_(a;-"-iE)D+(aj"+iE) 

_^(n+l)(/+l/2)^_^^n+l^)^^^^-n-l^) , (5.57) 

which is precisely ()3.74p . 

Next we evaluate (I5.57P at / = and n = M for integer M, and replace (7^*^^ with 
Tm/2 to find 

'-Tm/2{E) = D^{-E,0)D+{-E,0) ^ D{-E) . (5.58) 

The final equality follows from recalling that D-{E,0) and !)+(£', 0) are the even and 
odd subdeterminants for the spectral problem defined on the full real axis. This, up to 
a factor of z/2 arising from our normalisation of D{0), is precisely the original result 
conjectured in [1]. 

In fact, the quantum Wronksian relation ()5.56p was the initial source of the ODE/IM 
correspondence [1]. As a bridge between the ODE and IM worlds, it has some advantages 
over the 'TQ' approach we have stressed in the foregoing. In particular, it allows for a 
complete proof of the correspondence [2] . 
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5.5 Numerical techniques 

One of the byproducts of the ODE/IM correspondence was the reahsation that energy 
levels for Schrodinger problems can be calculated using nonlinear integral equations [1] 
of either TBA type or an alternative type known as the Destri de Vega or Kliimper- 
Batchelor-Pearce equations (or simply NLIEs). In appendix |D] we derive the TBA equa- 
tions and NLIEs associated with the spectral problems above, and explain how to obtain 
the spectrum by solving these equations iteratively. 

These methods of solving such eigenvalue problems appear to be new, though it- 
erative solution methods based on functional relations for spectral determinants had 
previously been employed by Voros [21], in work which was an important input to the 
initial observation of the ODE/IM correspondence in [1]. Numerically, the method is 
rather efficient - integral equations tend to be easier to solve than differential ones, and 
spectral determinants encode all the eigenvalues at once. 

5.6 The full dictionary 

All good correspondences need a dictionary, and the table below summarizes the mapping 
between objects seen by the integrable model and the differential equation. 



Integrable 




Schrodinger 


Model 




equation 


Spectral parameter 


<— > 


Energy 


Anisotropy 


<— > 


Degree of potential 


Twist parameter 


<— > 


Angular momentum 


(Fused) transfer 
matrices 




Lateral spectral problems defined at |a;|=cx) 


Q operators 




Radial spectral problems linking |.t|=oo and 
|xhO 


Truncation of the 
fusion hierarchy 


<— > 


Solutions on finite covers of C\{0} 



Armed with the dictionary, the horizontal axis of figure [2] can be annotated to indicate 
which integrable models correspond to the various values of M in the Bender-Boettcher 
problem. Thus for M = ^,1,|,2 and 3, the relevant integrable models are the A^=2 
SUSY point of the sine-Gordon model the free-fermion point, the Yang-Lee model, Z4 
parafermions and the 4-state Potts model respectively. It is amusing that the potential 
is related by the correspondence to the Yang-Lee model (or, strictly speaking, to the sine- 
Gordon model at the value of the coupling which allows for a reduction to Yang-Lee), 
thus returning by a very indirect route to the original thought of Bessis and Zinn-Justin. 
Note that while the various results discussed here are derived for M > 1, the 'Airy case' 
(M = 1/2, / = 0) has already been studied in [168], and has been shown to be consistent 
with integrable field theory predictions [1,169]. 

We finish this section with the table below, which records the notation used for the 
related objects appearing in the six vertex model, the continuum integrable model and 
the differential equation picture. 
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Six vertex 


Integrable 


Schrodinger 




TTl (\(\ pl 


equation 


u 


s 


E 






M 






I 




T 


C(2m) 


qQ{u,±<l)) 


Q± 
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6 Applications and generalisations 



The correspondence between ordinary differential equations and integrable models is in- 
triguing, but is it good for anything? One answer comes from the fact that the properties 
of ordinary differential equations and their solutions are rather better understood than 
those of the T and Q functions from integrable quantum field theory. The correspondence 
therefore allows a number of conjectured properties of these functions to be proven for 
the first time - analyticity, duality. . . [2, 153-155, 170-172]. 

On the other hand, ideas from the theory of integrable models have led to new insights 
into the spectral properties of a number of quantum-mechanical problems. This section 
begins with some examples of this aspect, starting with the reality properties in VT- 
symmetric quantum mechanics that were described back in section [2l First we sketch 
how the most general class of problems discussed there can be linked with integrable 
models. 



6.1 Inhomogeneous potentials 

A key feature of the proof of the TQ relation in section H] was the fact that the differential 
equation mapped to itself under the 'Sibuya' variable change x i— > io~^x, 
where u;= exp(7rz/(M-|-l)). In turn, this relied on the 'potential' V{x) = x'^^^ being 
invariant under V{x) ^ uj"'^ V {uj~^ x) . A natural generalisation, first observed in [7] 
and further explored in [11], is to add a term which is exactly anti-invariant under the 
same transformation. This leads to the consideration of the inhomogeneous potential 
V{x) = -|- ax^^~^ and the ordinary differential equation 



- dx"^ 



^{x) = E^{x), (6.1) 



where a is a free parameter. Again the radial spectral problem asks for those values of 
E admitting solutions which vanish as x ^ +oo, and behave as x'"*"^ as x ^ 0. With 
the inclusion of the additional term the subdominant large-x asymptotic, defining the 
fundamental Sibuya solution, becomes 



^ -M/2-a/2 , 



y{x,E,a,l) exp(^-|j^j. (6.2) 

For integer k we then generate a set of solutions to (|6.ip 

yfc(x,^,a,/) =a;'=/2+W2y(^-fc^^ ^2fc^^ ^(Af+i)fc^^ ^ a; = , (6.3) 

with each pair {yk,yk+i} forming a basis of solutions. Therefore, just as before, there 
must be a Stokes relation of the form 

C{E, a, l)yo{x, E, a, I) = y-i{x, E, a, I) + yi{x, E, a, I) , (6-4) 

with Stokes multiplier C{E, a, I) 

C{E,a,l) = W[y.,,yi]/W[yo,yi] = T^[y-i ,yi] . (6.5) 

Since uj^~^^ = —1, the only values of a to arise in the set of functions yk{x, E,a,l), 
related by (|6.4|) . are a itself, and —a. This makes it convenient to define 

CW(S) := C{E,±a,l) , y^^\x, E) := y{x, E,±a,l) (6.6) 
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after which taking (j6.4|) at positive and negative a imphes the pair of coupled equations 

=C^-(l+-)/2y(-)(c^X,t^2M^) ^^(l+a)/2y(-)(^-l^^^-2Af^^ (6.7) 

To remove the x-dependence we project onto solutions ip±{x,E,a,l) which are defined 
at the origin by 

^/;(x,a,/) = V+(2;,a,/) ~ + . . . , x^O (6.9) 

and ip-{x,a,l) = tplx, a, —l—l). The subscripts ± should not be confused with the 
bracketed superscripts (it), which are used to differentiate positive and negative values 
of a. We set 

D{E, a, I) = W[y, V'+J {E, a, /) . (6.10) 
Then, reasoning as before, the Stokes relations (16. 7[) and ()6.8p imply 

C'^+\E)D^+HE) = u;-(2i+l+a)/2^(-)(^2M^) ^^(2/+l+a)/2^(-)(^-2M^) (g^-^^) 
C(-) {E)D^-^ (E) = ^-(2'+l-«)/2^(+) (^2M ^) ^ ^(2i+l-a)/2^(+) ^^-2M ^) (g_^2) 

where D^^\E) = D(E, ita, /) are easily identified with the spectral determinants for the 
radial problems ()6.1|) taken at ita. These equations entwine the spectral problem at +a 
with that at —a, and reduce to two copies of the six-vertex TQ relation for a = 0. At 
general a [7] they match equations satisfied by the T- and Q-operators for the 3-state 
Perk-Schultz lattice model [173-175], a model with Uq{gl{2\l)) symmetry. 

6.2 PT-symmetry and reality proofs 

For maximal generality, we consider the three-parameter family of Hamiltonians 7iM,a,i 
and look for solutions of the corresponding Schrodinger equation 



.i.^_(,^)2M_^(.^)M-l^i([+l) 

dx'^ 



'ipkix) = Xki'kix) (6.13) 



that decay simultaneously at both ends of a contour C which can be taken to be the real 
axis if M<2, suitably distorted to pass below the origin if /(/ + 1) 7^ 0. For larger values 
of M the contour should be further distorted so as to remain in the pair of Stokes sectors 
which includes the real axis at M=l, as illustrated in figure [12] of section [H We shall 
prove 

Theorem: The spectrum of 7iM,a,i is 

• real for M > 1 and a < M + 1 + |2/ + 1| ; 

• positive for M > 1 and a < M + 1 - |2/ + 1| . 

This result includes the cases considered in section [2j 

First we remove the factors of i by setting <I>(x) = ip{x/i), so that the Schrodinger 
problem becomes 



dx^ x^ 



^kix) = -Xk ^kix) , $fc(x) G L\iC) , (6.14) 
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and has the same form as (|6.ip with E = —\k, though with different boundary conditions: 
to quahfy as an eigenfunction, <I> must decay as |x| — > oo along the contour iC. However, 
it is an easy generahsation of the discussion in section [5] that the function C{—X, a,l) 
defined in (j6.5p is the spectral determinant associated to the spectral problem (j6.14p . 
This identification will allow us to prove the theorem using techniques inspired by the 
Bethe ansatz. 

We start from equation (j6.1ip : 

and define the zeros of C^^\E) = C{E,a,l) to be the set {— A^}, and the zeros of 
D'^^\E) to be the set {e[^^}. Setting E = e[^^ in IKTE\\ leads back to a (coupled) set 
of Bethe ansatz equations for the zeros of D^^\ much as before. However we can also 
investigate the effect of setting E = —Xk- The left-hand side of (|6.15p is again zero, so 
the equation can be rearranged to read 

D'^-\uj^Xk)/D'^~\uj-^Xk) = . (6.16) 

For M > 1, WKB estimates show that the function D^^\E) has order less than one, so 
that it can be written as a simple product over its zeros, as in (I5.22p for the homogeneous 
potentials. Using this factorised form gives the following set of constraints on the A^'s: 

n = ' ^ = 1> 2, . . . . (6.17) 



n=l 



, E^n^ + W-2 X 



Notice what has been achieved here - the little-understood set of numbers {A^}, eigen- 
values of a non-Hermitian eigenvalue problem, are being constrained by the much better- 
understood set {En ^} . Since the original eigenproblem (j6.13p is invariant under / — > 
— 1— we can assume I > —1/2 without any loss of generality. Then each En ^ is an 
eigenvalue of an Hermitian operator TCM,-a,h ^-nd hence is real. Furthermore a Langer 
transformation [176] (see also [2,4]) shows that the En solve a generalised eigenproblem 
with an everywhere-positive 'potential', and so are all positive, for a < 1-I-2L This can 
be sharpened by considering the value of ■* (-^)I_b=o ' found by an easy generalisation 
of the result obtained in appendix iBl for the cases with a = 0. From 

^ \ 2M+2 "I" 2 y 

we see that D^^\E)\^^q first vanishes when a = M+21+2. Until this point is reached, 
no eigenvalue En ^ can have passed the origin, and all must be positive. (It might be 
worried that negative eigenvalues could appear from E = —oo , but this possibility can 
be ruled out by a consideration of the Langer-transformed version of the equation.) 

Taking the modulus^ of (16.171) , using the reality of the E^^ ^ , and writing the eigen- 
values of (|6.13p as Ajt = |Afc| exp{i6k), we have 



fr nEi^? + \Xk\' + 2Ei^\Xk\cos{,,^, , J ^ 
M V(4"^)2 + |Afc|2 + 24-)|Afc|cos(^ .^l • 
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For a < M + 2/ + 2 , all the En ^ are positive, and each single term in the product on 
the LHS of ()6.19p is either greater than, smaller than, or equal to one depending only 
on the relative values of the cosine terms in the numerator and denominator. These are 
independent of the index n. Therefore the only possibility to match the RHS is for each 
term in the product to be individually equal to one, which for 7^ requires 

cosij^ + Sk) = cosij^ - 6k) or sin(^) sin(,5fc) = . (6.20) 

Since M > 1, this latter condition implies 

5k = miT , m G Z (6-21) 

and this establishes the reality of the eigenvalues of (j6.13p for M>1 and a<M + 21 + 2 
or, relaxing the condition on a<M + 1 + |2Z+1|. 

The plots and discussion in section [2] indicate that most of the become complex 
as M falls below 1, at least for a = 0. We now see that this coincides with the point at 
which the order of D^~\E) is greater than 1, the factorised form of D^~\E) provided by 
Hadamard's theorem no longer has such a simple form, and the proof just given breaks 
down. 

The borderline case M = 1 is the simple harmonic oscillator, exactly solvable for all 
/ and a in terms of confluent hypergeometric functions. The case a = is discussed in 
detail in appendix IB.lj since a nonzero value of a can be absorbed in a simple shift of 
E, it is easily seen that the correctly-normalised solution for the general case is 



y{x, E, a, I) = ^x'+^e--'/^U{\{l + |) - / + |, x^) , (6.22) 

and that 



27r 

C{E, a, OIm-1 = V • (6-23) 



2 ' 4 j"- \2 

Thus the eigenvalues of (f6T3]) for M = 1 are Afc = 4fc - 2 - a ± (2/+1), A; = 1, 2, ... . 
They are all real for all real values of a and and are all positive for a<2 — |2/+1|. 

To discuss positivity at general values of M > 1, we can continue in M, a and I 
away from a point in the region {M=l, a<2 — |2/+1|}. So long as a remains less than 
M + 1 + |2/+1|, all eigenvalues will be confined to the real axis during this process, and 
the first passage of an eigenvalue from positive to negative values will be signalled by the 
presence of a zero in C(— A, a, Z) at A = 0. Using (j6.18p we have 

M+i\ atTT 27r 



C{E,a,l)\p_r,= {^Y'^—, f-^ ^, (6.24) 



2 + 2M+2 ) '' \1 2M+2 

and so the first zero arrives at -E = — A = when a = M + l — |2/+1|. Thus for all 
a<M + 1 — |2Z+1| , the spectrum is entirely positive, as claimed. 

Referring to the regions A, i?, C and D shown on figure [TSl the proof establishes 
reality for (a, V) & B (J C U D, and positivity for (a, I) G D. This doesn't mean that the 
spectrum immediately acquires complex elements when the region A is entered, though. 
In [12] the shape of the domain of unreality was investigated in more detail, for the 
special case M = 3 . (This value was chosen because the PT-symmetric problem is 
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Figure 18: The approximate 'phase diagram' at fixed M. The proof in the text shows 
that the spectrum is entirely real in regions B, C and D, and positive in region D. 

there equivalent to a radial problem, making the numerics a little easier to handle. This 
equivalence will be explained in the next subsection.) 

FigurefTOlrepeats the plot from [12] already shown in the prelude, this time superimposing 
the boundaries of regions A, B, C and D. The full domain of unreality (the interior of 
the cusped line), only touches the borders of region A at isolated points. There is also a 
small, approximately-triangular region inside A near p = within which the spectrum is 
entirely real and positive, despite being outside region D. 

In fact, the condition a<M -|- 1 -|- |2/ -|- 1|, which arose as a technical point in the 
proof, has a rather more physical explanation [12]. Along the lines a = M -|- 1 it (2/ -|- 1) 
the problem can be reformulated in terms of a "PT-symmetric version of supersymmetric 
quantum mechanics. On the supersymmetric lines the model has a protected zero-energy 
state which permits eigenvalues to become degenerate and pair off into complex conjugate 
pairs. For M = 3 this occurs at the points where the cusped line touches the borders of 
A on figure [T9l 

The idea of a protected zero-energy level also leads to a partial explanation for the 
pattern of cusps seen in figure [T9j In addition to the supersymmetric lines a = M+l it 
(2/-I-1), there is a protected zero-energy state along all line segments shown in figure [19) 
this follows from (16.240 . which shows that C{E, a, I) vanishes whenever a = {M+l){2n + 
1) lb (2/-I-1), n = 0, 1, . . . . For n > 1 these are precisely the lines on which cusps would 
be expected to be found, as follows from a consideration of the mergings of the relevant 
levels, and the fact - special to M = 3 - that the merging levels enjoy an E ^ —E 
symmetry. 

The cusped line marking the frontier of the zone of reality is in fact only the first 
of infinitely-many such lines - figure [20] shows some further lines across which pairs of 
energy levels become complex. 

Finally, we remark that Shin has extended the Bethe-ansatz-inspired ideas of the 
main proof given above to establish spectral reality for a somewhat wider class of non- 
Hermitian potentials [89,92]. 



51 



20 



10- 



a 



-lOJ 



-10 



10 



20 



Figure 19: The domain of unreality of 7is^a,p, where p = V3{21+1). Also shown are 
segments of the lines along which the problem has a protected zero-energy level. 



6.3 Curiosities at Af = 3 

The correspondence with Bethe ansatz systems can also reveal unexpected relationships 
between the spectral properties of apparently very different problems. One set of exam- 
ples, obtained in [11], connects quantum-mechanical problems via a third-order equation. 
Consider the radial problem studied in section [UTI at M = 3. For later convenience, we 
also swap / for the variable p := \/^{2l+l) already used in figure [T9l so that the differential 
equation is 



+ ax" + 



12 



^{x) = E<^{x) . 



(6.25) 



The spectrum is encoded in the spectral determinant D{E, a,l)\M=3- 

On the other hand, we can consider a third-ovdei ordinary differential equation: 



[D{g2-2)Digi-l)D{go)+x^]cPix) 



E(t>{x) 



(6.26) 



where D{g) = {-^ — Third-order equations of this type are generally associated with 
5C/(3)-related Bethe ansatz equations, as shown in [5,8], but for the particular 'potential' 
x^ that appears in (I6.26p . they collapse onto SU {2)-type equations [11]. With the choice 
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l+{a+V3p)/4: , gi = l+a/2 , 52 = l+(a-\/3p)/4 , 



(6.27) 



the collapse is onto the same set of equations as govern the problem (I6.25p . This leads 
to a rather remarkable spectral equivalence between pairs of second- and third- order 
ordinary differential equations. Furthermore, the third-order equation is symmetrical in 
{go^ 91t 92}, a feature which is completely hidden in the original second-order equation. 
By playing with this symmetry, some novel spectral equivalences can be established 
between different second-order radial problems, and also between these and certain lateral 
problems [11]. The parametrisation in terms of a and p allows these relations to be 
expressed in a compact way in terms of certain 2x2 matrices in the Weyl group of 
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Figure 20: An expanded view of figure [T9l showing the frontiers of the regions with zero, 
two, four and six complex levels. 

SU{3) [100]. With Q = ^ ) , let 

Radial(Q:) = Spect ( (j6.25p with radial boundary conditions ) 
Lateral(Q;) = Spect ( (j6.25p with lateral boundary conditions) 

and define matrices L , T and LT by 

-Kaf)' --K-".^^)- 

The reflections L and T together generate the Weyl group of SU{3), with their product 
LT being an anticlockwise rotation in the {p,a) plane by 27r/3. The three mappings 
are illustrated in figure Ell With this notation in place, three simply-stated spectral 
equivalences hold: 

(a) Radial(Q;) = Radial(Ta) 

(b) Lateral(Q:) = Lateral(LQ:) 

(c) Lateral(a) = Radial(LTQ;) . 

Of these, the second is trivial - it just encodes the I —1 — I symmetry of the 
lateral problem - but the other two are not. The third, (c), gives a direct insight 
into spectral reality for M = 3: it demonstrates that the PT-symmetric problem for 
this case is equivalent to an explicitly Hermitian problem. Notice that the domain of 
unreality for the PT-symmetric situation - shown in figure [19] above the wavy line - 
is rotated under (c) into a region of negative p, where I is sufficiently negative that 
the radial problem Radial(LTQ:) is also non-Hermitian and so no contradiction with 
unreality arises. (Interestingly, there are similar equivalences [110,177,178] for a subset 
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Figure 21: The three mappings involved in the spectral equivalences. 



of the non-Hermitian M = 2 problems, though the methods of proof in those papers are 
very different.) 

In the (a, p) coordinates adapted to M = 3 the lines with a protected zero energy 
level shown on figure [19] are a = 4 J it p/\^, J = 1, 3, 5, . . . , and they have an extra 
significance [12]: they correspond to points at which the model exhibits a PT-symmetric 
version of the quasi-exact solvability discussed in [179, 180], with an odd number J of 
'QES' eigenfunctions and eigenvalues which can be found algebraically. In fact, a further 
spectral equivalence from [11], related to a higher-order supersymmetry [11,181,182], can 
be used to show that if a lies on the line a = AJ + p/^/S, with J an integer, then the 
non-QES part of the PT-symmetric spectrum Lateral (o:) is the same as the full VT- 
symmetric spectrum Lateral(TQ;). (Likewise, if a lies on a = 4 J — p/\/3 then there is a 
mapping of the non-QES spectrum onto Lateral(LTLQ:) .) If a is inside the region A, so 
that some levels of Lateral(Q;) might be complex, then, looking at figure [2Tb. Ta must be 
outside this region. Thus the non-QES part of Lateral(Q;) is real, and so the levels which 
become complex as a moves into region A through QES values must lie in the solvable 
part of the spectrum [12]. This complements a similar, but as yet unproven, conjecture 
concerning quartic QES potentials [183]. 

6.4 Further generalisations 

The ODE/IM correspondence introduced here for SU{2) Bethe-ansatz systems has been 
extended to a number of other integrable models and ordinary differential equations. To 
finish, we present a very schematic summary of the integrable models and differential 
equations involved in correspondences to date. The cited references should be consulted 
for more details of the models and parameters involved in each case. 

• Six-vertex model, su{2) spin-^ XXZ quantum chains and the perturbed 
boundary sine-Gordon model [14]: 

As already described, the ground state of the six-vertex model is related to the 
spectral properties of the ordinary differential equation 




,2M 



+ 



/(/ + 1) 



E 'I'(x) = . 



(6.29) 
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This correspondence has been generahsed [14] to map the excited states of the 
integrable model to the following differential equations: 



.2M+2 



Zk) - E -^{x) = , (6.30) 



where the constants {zfc} satisfy 

L 

E 



^ Zk{zl + {M+2,){2M+l)zkZj + M{2M+l)z^) 



and 



{Zk - Zjf 



A 



Mzk 
4(M+1) 



+ A = 0. 



{21 + 1) 



AM' 



(6.31) 



(6.32) 



16(M+1) 

These equations were also obtained in [184]. Very recently, they and the simpler 
equation for the ground state have been interpreted as special cases of Langlands 
duality [185]. 

Perk-Schultz models and the perturbed hairpin model of boundary in- 
teraction [7,11,186]: 

Adding a inhomogeneous term to (|6.29p gives the ordinary differential equation 
discussed in section I6.lt 



dx"^ 



E 



■^{x) =0. 



Via a change of variables this is equivalently 



d 



-p^ + 2qtie^ + K^(e^^ + e" 



(6.33) 



(6.34) 



the form that was used in [186]. 

Spin-j su{2) quantum chains and the boundary parafermionic sinh- 
Gordon model [187-189]: 

The higher spin su{2) chains are found by a process known as fusion, and their 
ground state Bethe ansatz roots group into characteristic 'strings' in the imagi- 
nary direction of length 2j. This behaviour is captured by generalised eigenvalue 
problems of the form 



dx'^ x' 



or equivalently 



dx'^ 



■^{x) = 



*(x) = 



(6.35) 



(6.36) 



where Q = h + h . The formation of strings for this equation is discussed in [19]. 
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su{n) vertex models [5,6,8,19,170,187]: 

To encode Bethe ansatz systems related to su{n) with n > 2, it turns out to be 
necessary to go to higher-order ordinary differential equations: 



i-irDnig)-PKix,E)\^ix)=0, 

where g = {go, 91 ■ ■ ■ 9n-i} is a collection of 'twist' parameters 
Dn{g) = D{gn-i - (n-1)) . . . D{gi - 1) L»(go) , D{g) = 

and 



dx X 



(6.37) 

(6.38) 
(6.39) 



Pk{x,E) = {x''^"'' -Ef 

with h = 71 + 1, the dual Coxeter number of su{n), and the integer K gives the 
degree of fusion of the vertex model. 



Finite spin-j XXZ quantum chains at A = ibi [15]: 



dx'^ cosh^ X 
with a ={m + l)/{2j + 2). 

Coqblin-SchriefFer model [171]: 



-I— + hi 
dx 



d^ N{N + l) M{M + l) 2' 
H r-o 1- <7 



sinh x 



-I— + hn 

dx 



■^{x) = 0, 



(6.40) 



^(x) =0. 



(6.41) 



Notice that, after a variable change, this is a particular case of (j6.37p . 
Circular Brane [154]: 

-3-2 + ^^e^ + '^^ exp(eJ') ^{y) = . 



Paperclip models [155,172]: 

d^ n . o 1. 



dx"^ 



P 



x\2 



+ K^(l + e^ 



^{x) = 0. 



(6.42) 



(6.43) 



Note, setting x = y — logn and taking the limit n — > 00 reduces this to (|6.42p . the 
circular brane. The paperclip models and the perturbed boundary hairpins of [186] 
both exhibit the conformal hairpin boundary condition in their ultraviolet limits, 
but the perturbations are different. 

Finite spin-i XYZ quantum chain (ofF-critical deformation of the A= — ^ 
XXZ model) [17]: 

d r 92 

6^^^(x,g)=[-^ 



+ 9n(n + l)p{3x\q^) + c{q, n) 



(6.44) 



where P is the Weierstrass elliptic function, q = e*'^'^, r is the modular parameter, 

c((?, n) = -Mn + ^) ^\lo q3^ ' (6-45) 
and T?i is the elliptic theta function. 
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so{2n) vertex models [19]: 



(6.46) 



where g = {go,gi ■ ■ ■ gn-i} is a set of twist parameters as for su{n), 

= {n—l—gQ, n—l—gi, . . . , n—l—gn~i} , and the generahsed potential term Pk{x, E) 
is as in (j6.39p . but with h = 2n—2. 



so{2n + 1) vertex models [19]: 

where Pk{x, E) is as in (j6.39p with h = 2n— 1. 
sp(2n) vertex models [19]: 



Dn{g^)[^D^{g)-PK{x,E)(^] Pk{x,E) 



^ix,E,g)=0, (6.47) 



^(x,S,g)=0, (6.48) 



where Pk{x-, E) is as in (j6.39p with h = n+1. 
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7 Conclusions 



The main conclusion of this review can be stated very simply: it is that the T and Q 
operators that arise in certain integrable quantum field theories encode spectral data. 
This gives a novel perspective on the Bethe ansatz, and also a new way to treat spectral 
problems via the solution of nonlinear integral equations. As an application, we have 
shown how Bethe ansatz ideas led to a proof of a reality property in PT-symmetric 
quantum mechanics. Techniques from integr ability also shed light on the way that reality 
is lost [16]. The correspondence has proved useful in the converse direction as well 
[153-155,171]. 

There are many further problems to be explored, and here we list just a few. First, 
it will be interesting to find out how many other BA systems can be brought into the 
correspondence, beyond those listed in section 16.41 

Second, most of the correspondences established to date have concerned massless 
integrable lattice models, in a field theory limit where the number of sites, and of Bethe 
ansatz roots, tends to infinity. The first match to a finite lattice system has been reported 
in [15], linking the Poschl- Teller (Heun) equation in complex domain and spin--^^^ XXZ 
spin chains at special choices of A. It would be of interest to extend this result to other 
models. 

Third, while we have concentrated on describing the correspondence for the ground 
state eigenvalues of the T and Q operators, a spectral interpretation for the excited 
state eigenvalues ti(z^), ^2(1^) and so on has now been found for the massless SU{2)- 
related system in the continuum limit [14] . The complicated distributions of Bethe roots 
is reflected in the significantly more complicated nature of the Schrodinger problems. 
In particular the potentials are no longer real, even for the radial problems. Since all 
eigenvalues satisfy the same functional relations as for the ground state, it is expected 
that a correspondence for the excited state eigenvalues of all models should exist. 

Fourth, it would be of interest to extend the correspondence to encompass massive 
integrable models, a first step towards this goal having very recently been made in [17]. 

Finally, we should admit that our observations remain at a rather formal and math- 
ematical level. It is natural to ask whether there is a more physical explanation for the 
correspondence, but perhaps this question will be easier to answer once some of the other 
open problems have been resolved. 
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A The algebraic Bethe ansatz 



In this appendix, we fill in some details of the Bethe ansatz for the six-vertex model. 
The technique we describe is called the algebraic Bethe ansatz, an elegant formulation 
developed by the 'Leningrad school', which reveals the role of the Yang-Baxter equation 
in a particularly transparent way. Our presentation borrows heavily from [127], fixing 
one or two typos; other useful references are [190] and [191]. To keep the discussion 
self-contained, we begin by recalling some of the definitions from section [3l 

The local Boltzmann weights for the model are parametrised in terms of the spectral 
parameter v and the anisotropy r] as 



a(z^, r\) = sm{rj + iu) = W 



h{y, 7]) = sm{r] — iv) = W 



c(z^, rf) = sin(2r/) = W 
representing the weights W as in figure 
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Figure 22: The local Boltzmann weights. 
From these local weights, a 2^ x 2^ transfer matrix (z^) is constructed: 



{ft} 





{v)W 




{v)W 




{v)...W 








a2 


03 




OAT 



{v) (A.4) 



where a = (ai, 02 • ■ ■ cktv) and ex.' = (q'^, a'g . . . a'^) are multiindices. 

The task of the Bethe ansatz is to find eigenvectors ^ and eigenvalues t of T , so that 



(A.5) 



As mentioned in section [3l the first step is to make a well-informed guess as to the 
form of a putative eigenvector. The algebraic Bethe ansatz makes this guess with the 
help of an additional piece of machinery, the monodromy matrix T: 



{ft} 



a ^(32 


{v)W 


"/32"'^/33' 


{v)W 




{iy)...W 


/37v"'^ b 




"2 




03 


UN _ 



{u) (A.6) 
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The definition of Tab differs from that of T in the omission of the sum over one of the 
horizontal spins, (3i . Of course, T can immediately be recovered from Tab ! simply by 
performing this sum: 



(A.7) 



However, it turns out that important data is also hidden in the off-diagonal elements of 
T. As for the transfer matrix, illustrated in figure [8] of section [31 the definition T is most 
easily digested with the aid of a picture, figure [231 where in the second line we used the 
double lines as a shorthand for the entire collections of lines carrying the multiindices o: 
and a' . 



Tiv) 



ab 





a 




Figure 23: The monodromy matrix. 
For convenience, we write the components of T as 

T{u)ab = 

so that 



T{u)^^ T{u). 



Ail.) B{v) 



and, in pictures, 



C{u) 



T{v) = A{v) + D{u) , 



B{u) 



(A.8) 
(A.9) 




D{v) 




(A.IO) 



As matrices acting on 2^-dimensional vectors. A, B, C and D inherit an important 
property from the microscopic Boltzmann weights by virtue of the ice rule, which pre- 
serves the total fiux of arrows through each vertex and hence through each collection of 
vertices. If we define S to be the total number of up arrows minus the total number of 
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down arrows (the 'spin'), then A and D act on vectors with a definite value of S to give 
vectors with the same value, while B decreases S by two, and C increases it by two. 

We now come to a key property of the Boltzmann weights of the six-vertex model, a 
sufficient condition for integrability. There exists a collection of numbers R'^{v), making 
up the so-called R matrix and related to the Boltzmann weights via R{v) = W{v — irj), 
such that 



c 
. a . 



MW 



c J. 
6 



W 



5 ■ 
a c 
a . 



a c 





{v)Rt,{v-u') . (A.ll) 



This, a version of the famous Yang-Baxter relation, 'intertwines' the local Boltzmann 
weights at spectral parameters v and u' . It is best handled diagrammatically, representing 
R'^ as in figure [Ml The perhaps-awkward, but standard, placing of the indices reflects 



Kt(y) 




Figure 24: The R matrix. 

the interpretation of i? as a matrix acting on the tensor product Vi®V2, where Vi and 
V2 are the two-dimensional vector spaces, spanned by and <— , that are seen by the 
indices c and d respectively. The resulting pictorial representation of the intertwining 
relation is shown in figure [25l Apart from a shift in the spectral parameter, the entries 




Figure 25: Intertwining the local Boltzmann weights. 

of the R matrix are just the original Boltzmann weights, so in particular they respect 
the ice rule and the flux of arrows through the 'i?' vertex is conserved. Explicitly, the 
non-zero entries are 



i?ZjZj(i^) = iZ^^(z^) = v(YY = vfYj^ = a{v—iri,r]) = sm.{2rj+iv) 
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b{iy—iri,rj) = —sm{iu) 



c{i^—ii],i]) = sin(2r/) . (A. 12) 



It is also useful to write R'^{i^) as a 4 x 4 matrix, with the rows indexed by ab and the 
columns by cd , and both pairs of indices running — >— >, — ><— , < >, <— <— : 



( sin(27/+ii/) 

— sin(iz/) sin(2?7) 

sin(2r/) — sin(ii/) 

\ 



(A.13) 



sin(2?7+iz/) / 



This notation makes it easy to check the identity 



Rli{u)Rff{-u) = sm{2v+iu) sin(2r?-zz.) 



c xd 



(A.14) 



which will be used shortly. (In the parallel universe of integrable quantum field theory, 
this is related to a property called 'unitarity'.) 

Once the local intertwining relation is known, it is a simple matter to chain together 
N such equalities to find an analogous relation for the monodromy matrix. This is 
illustrated in figure [26l 




Figure 26: Commuting the monodromy matrices. 
In equations, with sums on repeated indices, 

Rt'i^^ - ^'){TXi.{^){T,-')M^') = {TXc{^'Wt')a'A^)R':^\i^ - . (A.15) 

It is usually convenient to leave the fact that each component of T is itself a matrix 
implicit, so that this becomes 

Rtiv - v')%^,{u)%,,.{u') = Tac{u')%.,.{v)Rt,{v - v') . (A.16) 
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Because of the way we have set things up, the ordering of matrices right to left in this 
equation corresponds to the ordering top to bottom of monodromy matrices in figure [26l 
Two important consequences now follow. First, multiplying both sides by i?^/^ {u'—u), 
summing on b and b' , taking traces on a, d and on a' ,d' , and then finally using the cyclic 
property of the trace and ()A.14p , it can be seen that the transfer matrices at different 
values of the spectral parameter commute: 



[TH,T(z.')]=0. 



(A.17) 



This establishes the claim made in section [3l and shows that all of the T(i/) can be 
diagonalised simultaneously, with z/-independent eigenvectors. 

The second consequence is more technical, but crucial for the algebraic Bethe ansatz 
approach to the construction of eigenvectors. The key idea is to treat 3(1/) as a creation 
operator, and for this it is necessary to know how to pass the other components of the 
monodromy matrix through it. Taking the relations implied by figure [261 over all possible 
values of the vector (a, a', 6, 6') of free indices gives a total of 16 quadratic relations. 
These can be rearranged to give, amongst others, the following exchange relations: 

: [B{u),B{u')]=0; (A.18) 
(^, ^, ^, H : D{u)B{v') = g{v-u')B{v')D{v) - h{v-v')B{u)D{v') ; (A.19) 

and (after swapping v and v') 

(^,^,^,^) : A{v)B{u') = g{u' -v)B{v')A{v) - h{v' -u)B{v)A{u') . (A.20) 

The 'structure constants' g and h are related to the components of the R matrix: 
, s R'Z-'Z^iv) a{h' — ir],r]) sin(2r/ + zz^) 



Rz 
Rl 



biv — if], 7]) 
c{ri) 



sin(zzv) 
sin(2?7) 



Rzl^{i^) b{v — irj,ri) sin{iu) 



(A.21) 
(A.22) 



One example should illustrate the derivation of these relations. Setting (a, a', b, b') = (— > 
,<—,<—,<—) in figure [26] and recording only the nonvanishing terms in the sums on left 
and righthand sides. 





+ v-v 




Translated into symbols, this is 

RZZ{i^ - y')D{v)B{u') + RZZ{t^ - u')B{u)D{u') 
which can be quickly rearranged to give (IA.19p . 



(A.23) 



B{v')D{u)RZZ{t^ -v') (A.24) 
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At last we can write down some eigenvectors. The first is straightforward: it is the 
false ferromagnetic ground state |r2) = | 1 1 . . . | ). In a vector notation, 

Since A and D cannot change the value of the spin S, and is the unique state with 
(maximal) spin N, it is an eigenvector of both A and D. In fact, 

A{i^)\n) = a^{u,7])\^}) , D{u)\n) = b^{u,7])\n) . (A.26) 

On the other hand C(i^)|il) = (there is no state with a larger value of S) , while i?(z^)|$7) 
is, in general, a new state. For the simplest example, at = 1, \^}) = (q) and A, B, C 
and D are all 2x2 matrices: 

Thus A{v)\VL) = a{v)\VL), D{v)\VL) = b{u)\Q), C{v)\^l) = 0, and B{v)\Vl) = c(i/)(°). 
Notice also that B[v)\Vt) would serve as a second eigenvector of T = ^4 + D. In fact, in 
this rather-trivial case it could equally have been B{v')\^), for (almost) any value of v' . 
To generalise this observation, at larger values of A^ one might search for eigenvectors of 
T of the form 

\^) = B{u{)B{u2)...B{unm. (A.29) 

This is the guess mentioned in section [3] as stage (i) of the Bethe ansatz. For A^ > 1 
this guess (or ansatz) does not always work - extra conditions need to be imposed on 
the numbers vi . . . z/„. To find these, we first compute the separate actions of A{v) and 
D{i') on our would-be eigenvector |^). Using the A—B exchange relation, (|A.20p . 

A{v)\^) = A{v)B{vi)B{v2)...B{vnm) 

= g{vi-u)B{ui)A{v)B{v2) . . . B{i^n)\^) 

- h{i^i-i^)B{i^)A{i^i)B{iy2) . ■ ■ B{i^n)\^) . (A.30) 

Continuing to pass A through the B operators, but only recording explicitly those terms 
resulting from the first term on the right-hand side of ()A.20p . yields 

n 

1=1 

n 

- h{ui-u)llg{u, - ui) a^(z., 7])B{u)B{u2) . . . B{un)\n) 

1=2 

+ further terms . (A.31) 

A nice way to reconstruct the further terms uses the observation that 1^*) depends on the 
parameters {i^i . . . Un} in a symmetrical manner, since, from (jA.lSp . the B{L'j) commute. 
This enables us to deduce that 

n 

A{um=A+\^)+Y,^t\^k) (A.32) 

k=l 
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where 



A+ = a^{u, 7]) Yl g{vj - u) , A+ = -a^{uk, v)K^k - oi^j - ^k) (A.33) 



and 



\^k)=B{v)\[B{v^m . 



(A.34) 



(It is instructive to check this directly, at least for n = 2.) 

The same steps can be repeated to find the action of the matrix D[u) on |^'), using 
(lAl9l) and ^K^. The result is: 



D[um = A-\^>) + Y.^'k\'^>' 



(A.35) 



k=l 



with 



A- = h\u, 7]) H g{u - Vj) , A- = -b''{vk.i^)h{v - v^) H ^(i^fc " ^o) ■ (^-36) 



The ansatz |^') will be an eigenvector of T = A + D if and only if all terms proportional 
to can be made to cancel between A(v)\^) and D{v)\^) . This requires A^ + A^ = 
for A; = 1 ... n, or 



sinh(2«r? - uj + Vk) b^iyk, v) 



k = 1 . . . n . 



(A.37) 



These are the Bethe ansatz equations (BAE) for the roots {^i.-.v^i}- It is important 
to reiterate that these equations do not have a unique solution, but rather a discrete 
set, matching the fact that the matrix T has many eigenvalues. For each self-consistent 
solution {vi} of the BAE, the corresponding eigenvector |^') of T has eigenvalue 



t{u) = A+ + A- 

n n 



(A.38) 



recovering the formula (j3.1ip quoted in the main text. 

The method can be straightforwardly extended to the more general situation of 
twisted boundary conditions [28,128]. The twist described earlier can be implemented 
by trading 



W 



aN 



and 



W 



{v) 



-i<f) 



w 



a J 



OtN 



(A.39) 
(A.40) 
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in the definition (|A.6p of the monodromy matrix, so that 'T{v)ab becomes 

_ (e-^^T^^ e^^T^^\ _ ( e-^^ A{v) e^^B{v)\ 

~ \e-^^r^^ e'^T^^J ~ \e-"^ C{u) e'^D{v) ) ' ^ > 

The more general transfer matrix 

T{u) = e-'*A{u) + e'^D{v) , (A.42) 

then has eigenvalues 

n n 

i=i i=i 
provided that the set of roots {z^i . . . z/„} solves the twisted Bethe ansatz equations 

jj^ smh(2i7? - Uj + Uk) {J^k,'n) 

Finally, note that the four-spin reversal invariance of the local Boltzmann weights 
implies the following property of the transfer matrix: 

T°'(z.,-0) = T|'(z.,</.) , (T=i,T=T) . (A.45) 

Indeed, from the definitions ([ATT]) . (IX2|) . ([O]) and (IXBl) we have T^' (i^) ab = T^' {T^)ab 
and therefore 

= T§V,<A) . (A.46) 

A first consequence of the equality ()A.45P is that T{u, (p) and T(z^, —(p) have the same set 
of eigenvalues. Furthermore, (|A.45|) also means that 

N 



T(zy, -</.) = C/T(z^,(/.) [7^1 , [/ = 0-^ (g) <T^ . . . O cj^ (A.47) 

where ax is in the vector notation a Pauli matrix. Hence = 1 and the set of eigen- 
vectors of T(i/, i;^) splits into [/-invariant singlets and doublets. Invoking the Perron- 
Frobenius theorem in the spin sector at (p = and continuity we conclude that the 
ground state is a singlet and its associated eigenvalue tQ{u, <p) depends on (p only through 
its square: 

to{i^,(P) = to{iy,-(P) . (A.48) 
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B ODE results 



This appendix collects some properties of a particular solution of the ordinary differential 
equation 



defined for M > 1 by the asymptotic 



y{x) = Ey{x) 



x-^^/^exp 



M+l 



as X tends to infinity in any closed subsector of S-i U 5o U 5i , where 

, . 2TTk 



arg^j; j 



2M+2 



< 



2M+2 



(B.l) 



(B.2) 



(B.3) 



This is the 'Sibuya solution', subdominant in Sq , relevant to the basic ODE/IM corre- 
spondence connecting the six-vertex model to an ordinary differential equation. 

B.l The special case of M = 1 

Setting M = 1 gives the simple harmonic oscillator, exactly solvable for all I in terms of 
the confluent hypergeometric functions M{a, b, z) and U{a, b, z). The correctly-normalised 
solution, subdominant in Sq, is 



y{x, E,l) = ^ x'+^e--V2f;( 1 + 3) _ |^ ^ + 3 ^ ^2^ 
V 2^ 



and has the large x asymptotic 
y{x,E,l) ~ 



/2i 



(B.4) 



(B.5) 



Notice the appearance of an extra factor of x^/^ compared with the formula ()B.2p . This 
has an interesting effect on the TQ relation [4j^. First, the modified asymptotic means 
that the most natural definition of the shifted solutions yk{x,E,l) is now 



yk{x,E,l) = w 



k/2~kE/2 



y{uj "x, uj^^E, I) , 



Jw/2 



(B.6) 



since this preserves the property Wq^i{E, I) = 1. The basic Stokes relation (j4.19p remains 

C{E, l)yo{x, E, I) = y_i(x, E, I) + yi{x, E, I) , (B.7) 

where C{E, I) = W-i^i{E, I). This must be projected onto a solution defined through its 
behaviour at the origin to eliminate x. Define 'ip{x,E,l) as in (15. 6p through the small-x 
asymptotic 

(B.8) 
(B.9) 



^{x,E,l) ~ x'+^ + 0(x'+'') 



and shifted solutions 



E, I) = a;'=/2+'=^/V(a;-^x, lo^'^E, I) . 



^^The argument here is a slight streamlining of that given in [4]. 
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The prefactor ensures that W[yk,ipk\{E,l) = W[y,ip]{uj'^''E,l) , while a consideration of 
the small-x behaviour of ipk shows that il)k{x^E,l) = uj^^^^^'^/'^'>^^^/'^il}[x^E^l). Hence 
W[yk,M{E,l) = u;Hi+m-kE/2 w[y,ij]{uj^''E, I). Taking the Wronskian of (fRTj) with ^ 
and setting D{E,l) = W[y,4^]{E,l), 

C{E, l)D{E, I) = a;-('+i/2)+^/2z)(cu~2^, /) + J'+^/^^-^/^D{u^E, I) . (B.IO) 

The modified asymptotic has resulted in a couple of extra E'-dependent factors. On the 
integrable models side of the correspondence, M = 1 is the free-fermion point and explicit 
constructions of the T and Q functions had led Bazhanov, Lukyanov and Zamolodchikov 
to precisely the same 'renormalisation' of the TQ relation [31], and so the correspondence 
survives. 

The functions C and D can also be calculated directly from the ordinary differential 
equation, of course. Taking from [192] the analytic continuation formula 



;7(a,6,ze2™) 
the Wronskian 

and, for 6 > 1, the \z\ 



(1 



r(i - b) 
r(i + o-&) 



M(a, 6, z) + e-2^^^"C/(a, 6, z) , (B.ll) 



W[U{a,h,z),M{a,b,z)] 
asymptotic 

r(6-i) 



m 

r(a) 



z-'e' 



U{a, b, z) 



Via) 



(B.12) 



(B.13) 



one finds 



and 



C(-E^>OIm=i 



D{E,l)\ 



2tt 



r(i + ™)r(i 



2l+l~E 



2r(/ + f 



M=l 



/2ir 



21+Z-E 



(B.14) 



(B.15) 



The zeros of D{E, I) are at -E = 2/ — 1 + 4k, k = 1,2, . . . , and reproduce the well-known 
eigenvalues of the radial simple harmonic oscillator. Zeros of C{E, I) correspond to there 
being normalisable wavefunctions on a contour which runs along the imaginary axis, save 
for a diversion to the right of the singularity at x = 0. A variable-change x — > x/i maps 
this to the real axis and results in the so-called PT-symmetric simple harmonic oscillator, 
discussed in [193] and [4]. Since this variable change also negates E, the eigenvalues of 
this problem, from (iBlIl) . are Ak - 2 ± (21+1), k = l,2,... . 



B.2 The values of C(0, /) and D{0, 1) 

For general values of M, the ordinary differential equation cannot be solved in closed 
form. However it does simplify at = 0, and this allows the value D{0, 1) to be found. 
Observe first that the function 



ip{x) = (^) 2*^+2 y ( [M+l'j M+l ^ MTT , E, I 



(B.16) 
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solves the Schrodinger equation 



+ X - CJX M+l + 



7(7 + 1) 



x^ 



ip{x) = A ip{x) , 



where 



M+l 



^ E , 7 = , A = . 

' ' M+l 2 ' 



(B.17) 



(B.18) 



At £'=0, a is zero and ()B.17p is the simple harmonic oscillator, which as just described 
is exactly solvable in terms of the confluent hypergeometric functions U{a, b, z) and 
M{a,b,z). The correctly-normalised solution, subdominant in the sector Sq, is 



\/2i 



(B.19) 



Reversing the variable changes, extracting the leading behaviour as x — > and comparing 
with ()5.13p . we find 



D{E, 01^=0 = ^r(l + ^)(2M+2) 
The value of C(0, Z) is easier to find: from (|5.14p at -E = 0, 



2/+1 

C(i?,0l^=o = 2cos(^^^ 



21+1 I 1 



2M+2^ 2 



(B.20) 



(B.21) 



B.3 The large-^ behaviour of D{E, I) 

Given that 'ip{x,E,l) is defined to have small-x behaviour ip{x,E,l) ~ the Wron- 

skian D{E,l) = W[y{x, E,l),^p{x, E,l)] can be evaluated for 5ReZ > —1/2 as 



D{E, I) = lim (2/+l)x' y{x, E, I) 



where 



y = 0. 



(B.22) 



(B.23) 



For this subsection the aim is to find the leading behaviour of D(E,l) as \E\ 00. As 
observed by Langer [176], the usual WKB approximation cannot be applied directly in 
this situation because of the singularity at x = 0. Instead, Langer suggested making the 
variable change x = , y{x) = e^/'^4>{z) . Setting A = / + 1/2 the equation becomes 



dz^ 



+ R{z,E,X) 







(B.24) 



where R{z,E,X) = e^^/^? - Ee^^ + and we set g = 1/(M+1). The general WKB 
approximation for (p is 



A 



R{z,E,X)^/^ 



exp 



y/R{u,E,X) du 



20 



+ 



B 



R(z,E,Xyl^ 



exp 



y/R{u,E,X)du 



(B.25) 

where zq is an arbitrary constant. Provided arg(— i?) < vr, this WKB expression is valid 
for all real z (and thus for x all the way down to zero) since 

i?-3/4(i2-i/4)" tends to 
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zero uniformly in z as — > +00 (see [158] for more discussion of the validity of the 
WKB approximation). The solution we aim to approximate is subdominant as x ^ +00, 
which requires A = 0. The value of B is fixed by the large-x asymptotic (IB.2p . It is 
convenient to set zq = +00, in which case the correctly-normalised approximate solution 
can be written for M > 1 {g < 1/2) as [4] 



■ exp 



^R{u,E,\) - e"/5 ]du-ge 



V2iR{z,E,Xy/^ 

Substituting into ([B?22]) gives the WKB approximation to the spectral determinant 



(B.26) 



lim 



' lim 

% 2— >— 00 



exp \z + 



du 



(B.27) 



To eliminate extraneous factors it is convenient to divide through by Z)^^^(0,/), giving 



Z)WKB(o,/) 



exp(/(i^,A)) 



where 



du . 



(B.28) 



(B.29) 



The remaining step is to analyse this expression as |£^| — > 00 with Qig{—E) < vr. Changing 
variables u^u + 1^ \u{-E) and setting = 1/(2-25-) = (M+1)/(2M) , 



i{E,x) = {-Erl_ 



00 
00 



du 



-EY 



M 



Performing the integral then gives the result quoted in the main text: 



lnD{E,l)r^^{-Er 



where 



ao 



r(- 



1 i_ 

2 2M 



|£'| ^ 00 , I arg(— £')! < tt , 



du 
dt{B.30) 

(B.31) 



^r(-^)r(^ + i) , /i = (M+i)/2M. 



(B.32) 



Notice that this result is independent of I, and applies equally to D^{E, I) = D{E, I) and 
(after analytic continuation) to D^{E,l) = D{E, —l—l). 
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B.4 The 'small- i?/large-/' asymptotic of In /). 

The Langer-transformed equation can also be used to extract the large-^ behaviour of 
the spectral determinant, a piece of data which is relevant for matching solutions of the 
quantum Wronskian equation. Once the initial differential equation for y{x,E,l) has 
been transformed into the equivalent Schrodinger problem ()B.24p . the role of the energy 
parameter is taken by — A^, and the WKB approximation is thus good in the large-Z limit, 
in addition to the large- -E limit investigated above. It is therefore possible to study the 
large-/ behaviour of D{E,l) with standard semiclassical techniques. We start with the 
WKB expression (Ib:28]) for the ratio D^^'^{E, l)/D^^^{Q, I), and expand the integrand 
of I{E, A) about = as 

V^eWg - Ee^^ + A^ - Ve^g + = - V ^ , — — , (B.33) 

so that 

oo 

I{E,\) = -Y,bn{X)E\ (B.34) 

n=l 

where the coefficients can be written in terms of standard Beta integrals: 



n!2" y^oo (e2e/9 + A2)"~^/^ 
^ ^i-2n+2ng ( Ting)T{-l/2 + n(l - g)) . (B.35) 

To ease comparison with[Cl we also give this result in integrable models notation, using 
the dictionary entries 

X = 2p/g , E = T\l-g)i-j s (B.36) 



to find 
where in the limit 1 



'^ 'ZZ^fn1} =-f:-n{p)s- (B.37) 



(r) \ 2n—2ng 
-] =a„p'-=»+'»9 , (B.38) 

»„ = rH)r(-i/2 + n(i-,)) ^, _ 

matching the results quoted in [153] and in ()C.2p . 
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C Quantum Wronskians and the Weiner-Hopf method 

In this appendix we show how the quantum Wronskian allows the power series expansions 
of the functions lnQ± to be pinned down uniquely, following Appendix A of [153]. A 
useful background reference for the Weiner-Hopf method is [194]. The proof relies on the 
following properties of Q{s,p): 

1. The function In Q^{s,p) has a formal power series expansion 

oo 

lnQ+(s,p) = -J]a„(p)s" . (C.l) 

n=l 

2. The coefficients o„(p) are meromorphic functions of p, analytic in the half-plane 
^e{2p) > —g, where g = (3"^ so q = e*'^^. 

3. As p ^ oo in the right half-plane, 

a„(p)~a„p^-2'^+'"^ (C.2) 

with 

a^ = -l^T{g)n\-g)T{l-g)\ (C.3) 

4. is related to Qj^ as 

Q.{s,p) = Q+{s,-p) = exp (^-f2an{-p)s^^ . (C.4) 

5. The quantum Wronskian condition holds: 

e^'''PQ+iqs,p)Q-iq-h,p) - e-^^'PQ+iq-\s,p)Q^iqs,p) =2ism{2TTp) . (C.5) 

These properties lead to a sequence of problems of Weiner-Hopf type which fix the coef- 
ficients a„(p) uniquely, as follows. 

First, substituting the power series (jC.ip and ()C.4p into the quantum Wronskian 
([05]) gives 



e'^^exp \^-Yliq''an{p) + g-"a„(-p))5" j 

-e-^-^'exp ^-^^(g-"a„(p) + g"a„(-p))s"^ =2isin(27rp) . (C.6) 

Expanding and equating coefficients leads to a set of relations of the form 

sm(TTng + 2'Kp)an{p) — sm{TTng — 2'iTp)an{—p) = Rn{p) , ?t. = 1, 2, . . . (C.7) 

where each function Rn{p) can be expressed in terms of ajt(ibp) with k = 1 . . . n—1. For 
example, 

Ri = (C.8) 

R2 = -{qai{p)+q-\i{-p)fe^^'Ps\n{2^p)/2. (C.9) 
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Let's consider a homogeneous case first, setting i2„ = 0. Then (|C.7p can be rewritten 

as 

r(l -ng + 2p) r(l -ng- 2p) 

Ting + 2p) ""^^^ = r(ng - 2p) ""^"^^ ^^•'"^ 

The left-hand side of this equation is analytic in the right half-plane 5Re(2p) > —g, and 
the right-hand side is analytic in the left half-plane 3^e(2p) < g. The common strip of 
analyticity —g < 3f?e(2p) < g ensures that both sides correspond to a single function 
fn{p) which is therefore entire: it is analytic on the whole complex plane. Furthermore, 
combining T{z + a)/T{z + 6) ~ z°'~'^ with the asymptotic ()C.2p shows that fn{p) = 
Odpp"^") in the right half-plane; and likewise in the left half-plane. For n > 1 Liouville's 
theorem then implies that fn{p) is identically zero (i.e. there is no zero mode) while for 
n = 1, is a constant, which is then fixed by (IC.3p . 

However, for n > 1 the RHS of ()C.7p is not zero. In addition the faster rate of decay of 
a„ (p) allows for a more effective rearrangement of (IC.7P . Start by rewriting this equation 
as 

sm{TTng + 27rp — ■KK)an{p) — sin^irng — 2-Kp — TTK)an{—p) = {—l)^Rn{p) (C.ll) 
for G Z, and then as 

r{l-ng + 2p + K) T{l-ng-2p + K) 

-an{p) = Sn{p) (C.12) 



r{ng + 2p- K) T{ng-2p- K) 



where 



Snip) = — r(l -ng + 2p + K) V{1 -ng-2p + K) Rn{p) . (CIS) 

vr 

This function has poles at the points 

2p = ±{K+l-ng + m) m = 0,1,2,... (C.14) 

coming from the gamma functions. Thus it has a strip of analyticity |3ftep| < [K+l — ng). 
To maximise the width of this strip, K should be as large as possible. At the same time, 
we must preserve the property that the first term on the LHS of ()C.12p decays in the right 
half-plane, and the second term decays on the left half-plane. The asymptotic behaviour 
of the first term in the right half-plane is 

SO the largest value we can take for K is n—2. The equation to be solved becomes 
T{n-l-ng + 2p) T{n-l-ng- 2p) 

■ (^n{p) - 777^ \ ^ an{-p) = Sn{p) (C.16) 



T{2-n + ng + 2p) T{2-n + ng-2p) 



where 



(-1)" 

Sn{p) = T{n-l-ng + 2p) T{n-l-ng- 2p) Rn{p) . (C.17) 

vr 

The final step is to find an 'analytic decomposition' of the right-hand side of (jC.16p into 
the sum of functions analytic and decaying in the left and right half-planes: 

Sn{p) = G+{p)-G^{p) . (C.18) 
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Assuming for the moment that this can be achieved, ()C.16p can be rewritten as 

r{n-l-ng + 2p) n t \ '^{n - I - ng - 2p) n i \ 

— — ^an(p) - G+{p) = — - G^{p) . (C.19) 

L [1 — n + ng + 2p) L [2 — n + ng — 2p) 

Now the same arguments as given above can be used to show that the left- and right-hand 
sides of (|C.19P are identicahy zero. Hence 

V{2-n + ng + 2p) 

an{p) = —TTT '^+(P> ■ (C-20) 

L [n — 1 — ng + 2p) 

It only remains to find the decomposition (IC.lSp . So long as Sn (p) has a non-empty strip 
of analyticity and remains bounded there, this can be achieved using Cauchy's theorem. 
Write 

where C is the contour aU bU cU d shown in figure [271 enclosing the point p and lying 
entirely within the strip of analyticity. 





d ^ 


\ Im 


P 










a 




c 




Re 




b 







Figure 27: The contour for the analytic decomposition of Snip)- 

Moving the segments b and d down and up to —ioo and +ioo respectively, their 
contributions to Sn{p) vanish, allowing us to write 

Sn{p) = G+{p)-G^{p) (C.22) 

where G-\-{p) is the contribution to ()C.21|) from a, and G^{p) is (minus) the contribu- 
tion from c, and these two functions do indeed decay in the right and left half-planes, 
respectively. Changing integration variable to t = p' /i and using (IClSh for Snip), (|C.20p 
becomes 



_ r(2-n-Fn5-F2p) r°° Tin-l-ng+2it) Tin-l-ng-2it) 

^niP) = 27r2 Tin-l-ng+2p) 7^ i^) ' 



(C.23) 
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and a„(p) has been determined uniquely in terms of the a^aip) with m < n. Since 
the first coefficient, ai{p), has aheady been fixed via the homogeneous equation and 
the asymptotic (1C.3|) . this completes the proof that the power series (IC.ip is uniquely 
determined by properties (i) ~ (v). Note that the question of the convergence of the series 
is at this stage moot. However, once the identification with the expansion of a spectral 
determinant has been made, this follows from standard properties of the solutions of 
ordinary differential equations. 
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D Derivation of the TBA equations and of the NLIE 



In this appendix we describe how to turn certain functional relations into nonlinear 
integral equations. We also indicate how the equations encode the effective central charges 
of the associated integrable models. 

D.l TBA equations from truncated fusion hierarchies 

Finite sets of functional equations such as the T-system (I3.63P - or equivalently (j5.45h 

- can be transformed into sets of nonlinear integral equations of the form known in the 
integrable model community as thermodynamic Bethe ansatz (TBA) equations [152]. 
Recall that the T-system (j3.63p is the following set of h—1 coupled functional equations 

t^"'\uj-'^E)t^"'\ujE) = 1+ Yl (t^^HE)) m = 1/2,1,..., (/i-l)/2 (D.l) 

i=i/2 

where uj = e27r«/ft,(M+i) q -g ^j^g incidence matrix of the Ah-i Dynkin diagram. In 
fact the following arguments apply equally for any simply-laced Dynkin diagram, and so 

- in the spirit of [195, 196] - we shall give a slightly more general treatment here than 
is strictly needed for the differential equation application. Also, because we consider the 
more general situation we shall use the integrable model notation T rather than t and s 
rather than E. We further streamline by setting r = h—1 and 

Tais)=&/^\E) , a = l...r (D.2) 

so that the T-system is 

r 

Ta{io-h)Ta{ios) = l + ll{n{s)f-\ a = l...r. (D.3) 

b=l 

In the more general cases, a and b label nodes on the Dynkin diagram described by Gab , 
and r is its rank. 

The first step towards the TBA is to introduce the Y- functions, defined by 

r 

Ya{s) = lln{sf^\ a = l...r. (D.4) 

6=1 

Clearly, impl les 

n{io-h)n{us) = 1 + n(s) (D.5) 

and a suitable product over b results in the set of equations 

r 

Ya{u;s)Ya{u;-h) = [](1 + ^(s))^'"' , (D.6) 

6=1 

called a Y-system. Equivalently, on setting s = e^/^ with = (M -|- l)/hM, we have 

h-l 

Ya{e + ii)Ya{e -ii) = J[{i + Y,ie)f"'' . (d.7) 

6=1 
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These coincide with the Y-systems of ALB. Zamolodchikov [195], which for the 
cases encode the finite size effects of certain integrable quantum field theories with Z/^ 
symmetry. (We shall comment on the relationship between these models and the contin- 
uum limit of the twisted six-vertex model at the relevant values of (p and r] at the end of 
this subsection.) 

Since in our cases T and Y are entire functions of s, they are 27ri/i = 27ri(/i+2)//i 
-periodic functions of 9. In fact, given the symmetry Ya = Yh-a, this periodicity also 
follows from the Y-system (]D.7P ; see [195,197,198] for more details. 

Next, the Y-system is transformed into a set of integral equations by introducing the 
' pseudoener gies ' 

ea{e) = lnYa{e) , (D.8) 

which as a result of ()D.7p satisfy 

r r 

ea{e + if) + ea{e -il)-Yl GabSbiO) = Gab Lbie) (D.9) 

6=1 6=1 

with 

La{e) = ln(l+e-^"(^)) . (D.IO) 

Equation (|D.9P has many solutions, and to pin things down some more information about 
the properties of the functions involved is required. We shall need the following. Since 
the functions t^"'/'^\s) are entire and hence regular at E = 0, the Ya{0) approach constant 
values as — > — oo . These constants solve the stationary (^-independent) version of ()D.7P : 
as will be justified shortly, the particular solution relevant here has all its components 
positive: 

sin2(^) 

On the other hand, as |s| — > oo with — vr + 5 < arg(s) < tt — 6 (where 6 is an arbitrarily- 
small positive number), the leading asymptotics of the ln(t^'")(s)), and hence also of 
the ln(ya(s)), are proportional to s'^lj. Thus, as ^e9 — > +oo in any strip |9m0| < 
7r(/i+2)//i — 6, the pseudoenergies behave as 

ea{e)r^maLe\ (D.12) 

with, in fact, exponentially-small corrections. A consideration of ()D.9P at large real values 
of 9, where the right-hand side vanishes, shows that the constants rUaL are proportional 
to the components of the Perron- Frobenius eigenvector of G. For the A^-i case of prime 
interest here, we can write these components as 

m^i = ^sin(^) , bo£R. (D.13) 

Finally, for the solution of (|D.9p that we require here - the ground state eigenvalue of 
the integrable model, or equivalently the quantum-mechanical spectral determinant ~ the 
zeros of the t^"'/'^\s) all lie on the negative s axis. This means that the Y-functions do 



^^Asymptotics of this nature can be proved in the ODE world using well-estabhshed techniques - see 
for example [20,89]; given the ODE/IM correspondence this is also now the most efficient approach for 
the integrable models. 
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not vanish on the real 9 axis, which when combined with ()D.12p justifies the choice of 
the solution (jP.lip of the stationary Y-system. It also implies that the functions 

fa{e)=ea{e)-maLe' (D.14) 

are bounded in the 'analyticity strip' |9m(0)| < tt/zJ^ and satisfy the relation 

r r 

fa{0 + il) + fa{e -q)-Yl Gabfm = ^ Gab ^{9) . (D.15) 

b=l 6=1 

Taking an (e-regularised) Fourier transform 

/CO 
d9f{9)e-''^'+^' (D.16) 
-oo 

of both sides of ([DlK]) . 

r r 
J^(2 6ab cosh(f ) - Gab) fb{k) = J2 Gab Hk) . (D.17) 
6=1 6=1 

Solving for fa{k) and transforming back to 0-space yields the equations 

1 . 

Sa{9) = maLe' - — Y, c^abi9 - 9')Lb{9') d9' (D.18) 
6=1 -^-"^ 

with 

r 

^ab{k) = -27r (2 5„,cosh(f ) - Gac)'^ G,^, . (D.19) 



c=l 

The inverse Fourier transforms of the kernels (pabik) can be found exactly in terms of 
elementary functions, and the results written as 

<Pab{9) = -i-^ In Sabie), (D.20) 

where 

Sab{0)= n W ' a,b = l...r. (D.21) 

The integer-valued index set Aab can be characterised via the Weyl group of the related 
Lie algebra, and 

{x} = (x-l)(x+l) , (x) = ■ (D.22) 

For Ah^i the explicit formula is 

a + b-l 

Sab{9) = H {x} , a,b = l...r. (D.23) 



\a-b\+l 
step 2 



^^The existence of such a strip surrounding the real axis is usually deduced numerically in finite-lattice 
models and extrapolated to the continuum limit. Via the ODE/IM correspondence this property has now 
been proven directly in the continuum: it follows from the entirety property of the T's and the fact that 
their zeros are exactly on the negative s-axis [11,89]. 
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For the other cases see, for example, [199,200]; for then' relation with the Fourier trans- 
forms (fDl9]) . see [195,196]. 

Nonlinear integral equations of the form (ID.ISP also arise in the framework of the 
thermodynamic Bethe ansatz [152, 201] method for finding the ground-state energies of 
massive relativistic integrable field theories in 1 -|- 1 dimensions, and hence they are often 
called TBA equations. In this picture L should be interpreted as the circumference of 
the infinite cylinder on which the theory is defined, and Sab{Q) are the two-particle S- 
matrix elements describing the factorised scattering of r particles of relative masses nia- 
Strictly speaking, the equations just derived correspond to an ultraviolet limit of these 
TBA equations, in which the overall mass scale is taken to zero. 

Once the ea{0) and hence the Ya{9) have been found via (jP.lSp . it remains to recover 
the T-functions. This can be done using much the same reasoning as led to the TBA 
equation, starting this time from ()D.5p . Dividing through by Ya{s), using ()D.4p and 
taking logs, 

r 

hiTa{e + il)+\nTa{e-il)-Y,Gab\rin{e)=La{e). (D.24) 

6=1 

From (jPTT]) and (|D^ . the funct ions InTa(^) — m^L e^/ (2 cos (vr/Zi)) are bounded in 
the same analyticity strip as the fa{&) above, and can likewise be found by a Fourier 
transformation. The final result is 

^''^'^^^^ = iZi/h) ''-hpjz ^'"^^ ~ ^'^^'^^'^ ^^'^^^ 

where the Fourier transform of the kernel ipab{Q) is 

i^ab{k) = -2tt (2<5,feCosh(f ) - Gab)'^ . (D.26) 

Again, the inverse transform can be performed, to find that tpabi^^) is defined as (pab{G)j 
save for the replacement of the blocks {x} by (x) in (|D.2ip . However, as explained in 
section IE. II below, for spectral applications one can often get away without calculating 
the Ta explicitly [1]. 

It is also interesting to extract the ground state energy of the spin chain, or of the 
corresponding quantum field theory. It might be tempting to use (I3.19P for this purpose, 
but one has to be careful as the continuum limit involves an infinite shift of the 0-like 
variable v in (j3.19p . This provides some intuition for the fact that the universal part 
of the continuum ground-state energy is to be found in the large-0 asymptotic of the 
fundamental transfer matrix, Ti{9) in current notations. 

The relevant asymptotic expansion has the form [30, 202] 

^ oo 



n=l 



where the numbers /Jn-i eigenvalues of the local conserved charges of the system 

on a cylinder of circumference L. The first of these, Ii, gives the ground state energy; 
and from the results of [30,202], 

lnT,(^) - ~ _ Asin{n/h) , _ _ _ 

^ ^ 2cos(7r//i) mi ^ ^ ^ 
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This asymptotic also follows from ()D.24p . Let q'a'^ be an orthonormal set of eigenvectors 
of Gab, where a runs over the exponents of G, so that 



Y^Gatq^^^ =2zos{Tia/h)q\^^ , Y.q\^) q^^^ = 5ab ■ (D.29) 



6=1 

Then 



^ cosh(7r/c/n) — cos[iTa/h) 

and is easily seen to have poles at A: = (ztcr + 2hn)i, n G Z. For large 6 the inverse 
Fourier transform yielding ipab{G) can be expanded over the residues of these poles, the 
leading term coming from the pole nearest to the real axis: 

MO) = — / dk^{k)e'^' 



27r 



e^^ + ... . (D.31) 



sin(^//i) 

For the A,, case relevant for comparison with the expansion (|D.27p , the normalised eigen- 
vector has components q^a^ = sin(7ra//i). Substituting ()D.3ip into the integral for- 
mula (1D.25P and swapping sums with integrals gives an asymptotic expansion which 
matches ()D.27p . and begins 

lnri(0) - ^ '^]^,^, e' = -y r sin(7r6//i) 6^'+'' ^{9') dO' + . . . . (D.32) 
2cos(7r//i) '^t^iJ-oo 

Comparing leading terms, 



— / mae'La{e)de (D.33) 

a=l -^-"^ 



where La{d) = ln(l + e^'^^'^'i) and £{6) solves (jDT8]l . Recalling that Cefi was defined via 

i^o(L) = F(L) = -^, (D.34) 
we see that ()D.33P is equivalent to 

dOmaLe'Laie). (D.35) 



r 

6 ^ 

CcS 



Even though £{6) cannot be found in closed form, it turns out that the integral in ()D.35P 
can be calculated exactly as a sum of Rogers dilogarithm functions 



1 r^' 

^x) = -- dy 



Iny ^ ln(l - y) 



(D.36) 
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(For more about sum rules for Rogers dilogarithm see [203].) For our model, the final 
result is 



Ceff 



a=l ^ 

h-1 2M - 1 



h + 2 2M + 2 




(D.37) 



which is exactly half the central charge for the Z/^ parafermion models [204]. How does 
this compare with the continuum limit of the twisted six-vertex model? As mentioned 
before the untwisted six- vertex model is equivalent, in its continuum limit, to the theory 
of a free compactified boson with central charge c^g = 1. For the twisted variant of the 
six-vertex model we have instead 



Inserting the relations 



we find 



TT M , TT /T^ X 

^=2mTI' '^ = ^mT^ 

rev _ 2M-1 

which is the result quoted in ()D.37p . The factor of two discrepancy with the parafermionic 
central charge is, for h odd, simply due to a double counting that can be avoided via a 
folding procedure that identifies conjugate nodes j ^ h — j in the related Afi_i Dynkin 
diagram. A slightly more subtle folding phenomena is at work in the h even case too [1]. 

D.2 NLIEs from Bethe ansatz systems 

The TBA method relies on the truncation of the fusion hierarchy, and so can only work 
at rational values of ry/vr or M. In addition the resulting equations depend in a com- 
plicated way on the arithmetic properties of these numbers. However, there is another 
approach which works more generally. Employing ideas developed in [5,27,28,31,205], 
the Bethe ansatz equations associated with a TQ relation of the type discussed here can 
be transformed into a single nonlinear integral equation. Starting with the TQ relation 
()5.14p . set E = and invoke the entirety of C{E, I) and D{E, I) to obtain 

-""^#^^ + 1 = ' ^ = 1'2,... . (D.41) 

Replacing D{E,l) with its Hadamard factorisation (j5.22p . valid for M > 1, the Bethe 
ansatz equations for the eigenvalues of D{E,l) are recovered: 



oo 



n 



Ej — oj'^Ek 



Ej - u-^Eu 



J — LU i^k \ _ 21 + 1 



uj''+^ , k=l... . (D.42) 



Setting 



"^^^-'^ D{u-^E,l)-'' 1\.[e,-u:-^e) ' 



81 



the Bethe ansatz equations are rephrased as conditions on the function a{E): 



a{Ek) + l = Q , k = l,2,... . (D.44) 

In fact, it fohows from the TQ relation that the zeros of a{E) + 1 are precisely the zeros 
of D(E), the Bethe roots Ef^., together with the zeros oiT{E). For simplicity, we consider 
the ground state situation, for which all the E^ are real and positive, and all the zeros 
of T{E) are real and negative. We begin by taking the logarithm of (1D.43P 

In a{E) = ^^^^ + f; F{E/E,) , (D.45) 
k=i 

where F{E) = ln[(l — oj'^E)/{l — uj^'^E)] . The logarithmic derivative 9£;ln(l + a{E)) 
has a simple pole at each eigenvalue E^. Applying Cauchy's theorem, the infinite sum 
in ()D.45P can be replaced by a contour integral, the contour encircling the points Ej. in 
an anticlockwise direction, while avoiding all other singularities of ln(l + a{E)). Given 
our assumptions about the locations of the zeros of D{E) and T{E), a suitable contour 
C runs from +00 to above the real axis, encircles the origin then returns to 00 below 
the real axis. With this in mind ()D.45P becomes 

lna{E) = ^^^^ + ^ ^ F{E/E') Oe' ln(l + a{E')) . (D.46) 

We change variables via E = exp(2M9/{M+l)) and define (with a mild abuse of nota- 
tion) In a{6) = lna(e^*^^/*^*^^^^) . Integrating by parts we have 



lna{e) = ^^l^-/ de'de>R{e-e')ln{l + a{9')) 
M+i Jc^ 

+ / dO' de'R{e - e') ln(l + a{e')) , (D.47) 

JC2 



where 



i?(0) = Aa^i7(e2A'/e/(M+i)) ^ p_48) 
27r 

The new integration contours Ci and C2 run from —00 to 00 just below and just above 
the real axis respectively. Using the property [a(^)]* = a{6*)~^, which follows from the 
Bethe ansatz equations, we rewrite the integrals in terms of integrations along the real 
axis: 



/oo 
d9'R{e -0')\na{e') 
-00 



m(2/+l) 
M+1 



/oo 
de'R{e-e')Qm ln(l + a{9'-i0)) . (D.49) 
-00 



It is then a simple matter to solve this equation using Fourier transforms. Using the 
notation f{k) and J^[f{9)]{k)) introduced earlier, the Fourier-transformed equation is 

(l--R(fc)).F[lna](A;) = -^^^^&i^5(A;)-2ifl(A;)9m.F[ln(l + a)](A:) . (D.50) 
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We apply (1 — R{k)) to both sides, then take the inverse Fourier transform to find 

/oo 
de' ip{9 - {9' - iO))9m ln(l + a{e' - iO)) , (D.51) 
-oo 

where 

ip{e) = T'^ [ (1 - R(k))-^R(k)m , (D.52) 

and mL is a real constant which arises from a zero mode and can be traced to the pole 
in (1 — R{k))~^ at k = i. It should be fixed by a consideration of the large-0 asymptotic 
of In a{9). From the definition of R{9) and the relations 

„ , sinh (t6' + zvrr 2cTsin27rT ,^ . 

iOe In — - — = — (D.53) 

smh a9 — ittt cosh 2a9 — cos 2itt 

d9 2asin27rr ^ sinh(l-2r)g 

27r cosh 2a9 — cos 27rr sinh 

2a 



oo 



we obtain a compact expression for the kernel '^{9): 

r dk sinh ^^(^'^'^ 

^(^) = / ^ e-^^ ' , . (D.55) 

' 7 27r 2sinh^cosh^ ^ ^ 

Finally we rewrite the equation in terms of the integration contours Ci and C2 ■ 

lna{9) = iTT{l + \) - imLe^ 



+ 



[ d9' ip{e-9')ln{l + a{9')) - [ d9' (p{e-9')ln{l + a{e')~^) . (D.56) 

JCi JC2 



This equation holds for all 9 within the strip |9m-0| < min(7r, vr/M). A little extra care 
is required for larger Qm9 as the kernel ip{9) has poles at 9 = ibzvr and itivr/M. If 
M > 1, the correct analytic continuation for positive values of Qm9 is given by the 
second determination [206] 

\na{9) = -i{l-e-'^/^)mLe^ 

+ [ d9' ^n{9-9')\n{l + a{9'))- f d9' ^n{9-9')\n{l + a{9')~^) , (D.57) 

JCi JC2 

where 

2icos(.5^)sinh(^- ^) 
^''^^ ^(cosh(2^-^)-cos(i^)) ^ ^ 

This equation is valid for vr/M < 9m < vr, when M > 1. A similar continuation can 
be performed for M < 1. For a detailed discussion of the effect of this continuation on 
eigenvalue asymptotics, see [16]. 

The NLIE (|D.56P first arose in [27] as the continuum limit of an equation describing 
the finite size effects of the six-vertex model. It also appears in relativistic integrable 
scattering field theory in 1 + 1 dimensions, where it describes the finite size effects of 
the ultraviolet limit of the massive sine-Gordon model [205]. In the latter situation the 
kernel (p{9) is related to the scalar factor of the sine-Gordon soliton-soliton scattering 
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amplitude. In the field theory context, the effective central charge is given in terms of 
In a{6) according to 



3imL 



Ceflf 



/ de e'^ \og{l + a{e')) - [ de \og{l + a-\e')) 



uci JC2 
This integral can be evaluated exactly, with the result 



(D.59) 



Ceflf 



6(/ + ^)^ 
M + 1 



(D.60) 



which, given the relations = 1/(M + 1) and p = {21 + 1)/(4M + 4), agrees perfectly 
with 



Ceff = c - 24Ap 



1-24^ 
/32 • 



(D.61) 



The above NLIE has been derived for the ground state, for which the Bcthe roots 
{Ek} are all real and positive. However, nonlinear integral equations of both TBA and 
NLIE types can be found even when some of the Bethe roots lie in the complex plane. 
See, for example, [31,202,206-209]. 
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E Calculating the spectrum of an ODE 



We briefly describe how to calculate the spectrum of 



E 



V'(x) = 



(E.l) 



using the nonlinear integral equations derived in appendix iDl 



E.l TBA approach 

When M is an integer and + =0 the spectrum encoded in the spectral determinant 
D(E, I) can be determined from the fusion relations ()5.45p using TBA equations. Recall 
the following relation obtained in section [3?71 

'-Tm/2{E) = '-C^''\v-'E) = D{-E) ^ D+{-E)D^{-E) . (E.2) 

From this we see that the real positive zeros of D(E) are precisely the real negative zeros 
oiTM/2{E) or C^^'^\E). We set 

t^'^\e + m^)=Tm{-E) , h = 2M. (E.3) 

Then from 

i + Y2n,{e + q) = t^"'\e + i7r^)t^^^{e + i7r^) , (e.4) 

we see that the eigenvalues {Ei^} are those zeros of 1 + Ym{0) that lie along the line 
3=7716* = 7r/2. 

The constants niaL = ^ sin ^ appearing in the TBA equations (ID.lSh must be tuned 
to match the asymptotic properties of the spectral determinants. The correct result is 
found by setting bo = 2cos(7r/2M) oq where oq was defined in ()5.20p . 

A very simple iteration scheme can be used to solve the TBA equations (jP.lSP nu- 
merically. Starting from 

ei"+i)(e) = a4-)(e) + (l-a) 

with e''a\e) = ruaLe^ and < a < 1, dES]) can be iterated until the desired accuracy 
(~ 10^^^) is reached. Empirically the value a = 0.5 gives a very good rate of convergence 
usually in less than one hundred iterations. Once the functions ea(^) are known numer- 
ically along the real axis, the TBA equations (jP.lSP provide an integral representation 
which can be used to reconstruct the functions ea{0) and Yai9) = exp(ea(6')) everywhere 
in the complex plane. The only point to watch is that the singularities of the kernels 
4>ab^ at ib77r//i and beyond necessitate the introduction of extra terms when the e^'s are 
continued beyond the strip — 7r//i < "^mO < n/h. 



TTlaLe 



-y 



MO 



(E.5) 



E.2 NLIE approach 

As mentioned above, the single nonlinear integral equation derived in appendix ID. 21 
encodes the zeros of both the spectral determinants C and D. The zeros of 1 -|- a{E) 
on the positive real axis are the zeros of D(E,l), while those on the negative real axis 
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are the zeros of C{E,l). Under the variable change E = exp(2M0/(M+l)) the positive 
real axis of the complex S-plane becomes the real-0 axis. Therefore the eigenvalues 
encoded in D{E, I) can be found by searching along the real axis for the zeros of l + a{9). 
Similarly the eigenvalues of C{E, I) can be found as zeros of 1 + a{9) on the line 9m 9 = 
7r(M+l)/2M. Before we can solve the NLIE (|D.56p we must fix the constants mL. The 
large- asymptotics of a{E) are 



log aiE) 



-liboil-e-'^/^^){EY 

iMEr 



2tt 



2tt 



M+T < arg(^) < 27r - ^ 
-^<arg(i?)<^ 



AI+1 

-2tt 



(E.6) 



M- 



- < arg(E) < 



M+l 



where bo = 2 cos{j^)ao as before. Above, by {E)f^ we imply e^^'^^^^^^\E\'-'- . Thus the first 
and third asymptotics coincide, as indeed they must since a is a single-valued function 
of E. Given (|5.27|) . we set mL = ^bQV~^, though the factor of v~f^ is arbitrary and is 
purely to match the conventions of [4,31]. 

The NLIE can also be solved by iteration and the function a{9) constructed for all 
complex 9 using the integral representations (|D.56P or (|D.57p as appropriate. Conse- 
quently the eigenvalues of the spectral determinants D{E) and C{E) can be calculated 
for all M > and I. This method is more generally applicable than the TBA approach 
simply because the NLIE depends on the parameters M and / in a continuous manner. 
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